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ABSTRACT 



Context. The Crab nebula and pulsar have been widely used as a calibration source for X-ray instruments. The in-flight effective area calibration 

of the Reflection Grating Spectrometers (RGS) of XMM-Newton depend upon the availability of reliable calibration sources. 

Aims. We investigate how the absolute effective area calibration of RGS can be obtained using Crab as a standard candle. 

Methods. We have analysed RGS observations of the Crab using different instrument configurations and spatial offsets, and made use of 

previous determinations of the continuum spectrum of the nebula plus pulsar. Due to the high spectral resolution of the RGS, we resolve the 

main absorption edges and detect the strong ls-2p absorption lines of neutral oxygen. 

Results. We get an excellent fit to the Crab spectrum using this fixed continuum and the absorption spectrum determined by RGS. We get 
accurate column densities for the neutral atoms of H, N, O, Ne, Mg, and Fe, as well as a clear detection of Fe ii and firm upper limits for other 
ions. Our data are in good agreement with earlier optical and UV spectroscopic measurements of some of these ions. We find solar abundances 
for N and O, while Ne is overabundant by a factor of 1.7 and Fe is underabundant by a factor of 0.8. We confirm that there is less dust in the 
line of sight compared to the prediction based on the absorption column. Our spectra suggest a more prominent role of ferric iron in the dust 
compared to ferrous iron. 

Conclusions. Our high-resolution observations confirm that Crab can be used as an X-ray calibration source. RGS spectra have determined 
the absorption spectrum towards Crab with unprecedented detail. 
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1. Introduction 

The Crab nebula and pulsar have been widely used as a calibra- 
tion source for X-ray instruments, starting already in the early 
days of X-ray astronomy. The main advantage of Crab is that it 
is one of the brightest steady X-ray sources on the sky and that 
it has a simple power-law spectrum. 

Later, more complexities of the spectrum have be- 
come clear. For instance, recent imagin g data (XMM- 



Newton: IWillingale et alJ 1200 It Chandra: IWeisskopf et al 



20001: lMoriet al.ll2004 show a complex spectral and spatial 



structure with the components such as pulsar, torus, jet, dust- 
scattering halo, etc. Also, the spectral shape of the pulsar and 
nebula are dif ferent, and both sh ow spectral breaks at higher 
energies (e.g.. lKuiper et al ] |200ll) . 

A complicating factor is that the interstellar absorption be- 
comes very high at low energies. For many instruments, the 
uncertainties in the low-energy calibration of the instruments is 
formally compensated for by adjusting the interstellar absorp- 
tion column density to get a best fit. This procedure, however. 
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introduces high systematic errors both in the derived eff^ective 
area at low energies, as well as in t he derived in t erstell ar ab- 
sorption. This issue is illustrated by iKirsch et al.l (120051) . who 
compare a large number of current and past X-ray instruments 
that have observed Crab. 

There is still an urgent need for accurate calibration of 
X-ray instruments in space. First, an accurate knowledge of 
the relative effective area is needed for sources with com- 
plex intrinsic continua, for example active galactic nuclei 
with warm absorbers, relativistically broadened iron or oxygen 
lines, soft spectral components and hard reflection components. 
Furthermore, an accurate knowledge of the absolute effective 
area is important when X-ray data are to be matched to opti- 
cal/UV data (for example for isolated neutron stars), or where 
the absolute flux has cosmological implications. A good exam- 
ple of the last issue are clusters of galaxies; their absolute flux 
determines their absolute hot gas mass as well as the inferred 
baryonic fraction. 

The problem with most previous and current X-ray instru- 
ments that have observed Crab is their relatively poor spectral 
resolution at the low energies where interstellar absorption is 
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important. The situation has improved considerably by the ad- 
vent of high-resolution grating spectrometers. With high spec- 
tral resolution, it is possible to resolve the absorption edges 
and absorption lines, thereby constraining the interstellar ab- 
sorption. With the Chandra Low Energy Transmission Grating 
Spectrometer (LETGS) it has been possible t o obtain a high 
resolution spectrum of the Crab pulsar ( Ten nant et al ]|2001), 
but the nebula containing most of the X-ray flux is too ex- 
tended for the LETGS to obtain a high resolution spectrum. On 
the other hand, the Refle ction Grating S pectrometers (RGS) of 
XMM-Newton dden Herder et al.ll200ll) are well-suited to ob- 
serve the full Crab nebula. 



From an astrophysical point of view, these Crab spectra of- 
fer us a unique opportunity to study the interstellar medium. 
For instance, the determination of the interstellar abundances 
from optical and UV spectroscopy is not straightforward, as 
only the gaseous components yield absorption lines. Moreover, 
these absorption lines are often heavily saturated, thereby mak- 
ing an accurate determination of ionic column densities diffi- 
cult. Furthermore, depletion of free atoms onto dust compli- 
cates the situation considerably. In X-ray spectra essentially 
all atoms, whether in gaseous or molecular form, contribute 
to the X-ray opacity, and the correction factors for differences 
in dust and gaseous opacity are not very high. Moreover, high- 
resolution spectra near the dominant absorption edges can help 
to disentangle the gaseous from the dusty or molecular phases, 
and to give clues on the chemical composition of the dust. 



In this paper we analyse a set of observations of the Crab 
nebula, obtained with the RGS of XMM-Newton. These data 
were taken for the purpose of effective area calibration. The 
paper is laid out as follows. In Sect. |2] we discuss the intrin- 
sic (unabsorbed) continuum spectrum of Crab as derived from 
other instruments. In Sect. [3]we describe the basic data extrac- 
tion and non-standard instrumental corrections that we need to 
apply. Sect.[3]then presents the results of our spectral analysis. 
In Sect, m we discuss our results on the interstellar absorption 
and the composition of the interstellar medium towards Crab. 
In Sect. |5]we discuss the implications for the calibration of the 
RGS and give a breakdown of the systematic uncertainties in- 
volved in this calibration. Our conclusions are given in Sect.|6] 



This paper is the fi rst of a series of three papers. In paper II 
dKaastra et al. I l2008bl) we present a new analysis of the X-ray 
model spectra of the white dwarfs Sirius B and HZ 43. We show 
in that paper how these models can be used to calibrate the long 
wavelength absolute effective area of the LETG S of Chandra 



with a n accuracy of a few percent. In paper III dKaastra et al 



2008ah we combine these results based on white dwarfs with 
our present results on Crab using blazar spectra taken simulta- 
neously with the RGS and LETGS. This leads to a refinement 
of the calibration presented in the present paper, and allows to 
determine the absolute flux of Crab within a few percent uncer- 
tainty. 
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Fig. 1. Measurements of the Crab nebula by iToor & Seward 
( 1974i) . Circles: proportional counter; stars: scintillation 
counter; filled symbols: data from Fig. 3 of Toor & Seward; 
open symbols: same data but corrected for Galactic absorption. 
Thick solid line: power-law model of Toor & Seward; thin solid 
line: our power-law fit to the same data; dashed line: power-law 
fit to absorption corrected data. 



2. The intrinsic spectrum of the Crab nebula and 
pulsar 

2.1. The Crab nebula as a calibration source 

have studied the Crab nebula as a potential calibration source 
for X-ray astronomy. They measured Crab in 1970 using a 
rocket experiment with proportional counters between 1.5- 
10 keV and scintillation counters between 5-70 keV. The mea- 
sured spectrum (nebula + pulsar) has a power-law with photon 
indexF = 2.10 + 0.03 and normalisation at 1 keV of9.7 + 1.0 
(in units of 1 0"* photons m"^ s^ ' ke V" ' , used throughout this pa- 
per for the Crab normalisation). Apparently no correction for 
Galactic absorption at the lowest energies was made. 

If we fit the data as presented by Toor & Seward, we find 
instead r = 2. 13 ±0.05 and = 10.3±0.9, i.e. a slightly steeper 
power- law. Correcting their data for Ga lactic absorption using 
simply iMorrison & McCammon dl983h absorption cross sec- 



tions and A^h = 3.45 x 10^^ m ^ dSchattenburg & Canizares 



.1986), we obtain a 20 % higher flux at 1.7 keV, the low- 
est energy measured by Toor & Seward. This results in F = 
2.18 + 0.05 and = 11.5 + 1.0, i.e. an increase in F of 0.05. 

Fig. [T] shows a comparison between these fits and the raw 
and absorption-corrected data of Toor & Seward. It is evident 
that at any energy the uncertainty in the spectrum is at least 
10 %. Also note the discrepancy between both instruments near 
10 keV: the scintillation counter measures about 13 % less flux 
in the common 5-1 1 keV band. 

From the above discussion it is clear that we need to pay 
more attention to the way the spectral parameters of Crab are 
derived and to their statistical and systematic uncertainties. As 
it is known that the pulsar and the nebula have a diff'erent spec- 
trum, we also have to derive both model spectra separately. 
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even though our RGS observations have insufficient spatial and 
time resolution to distinguish both. 



Since the early work of iToor & SewardI (119741) . many mea- 
surements of the Crab spectrum using differe n t instr uments 
have been made, see for instance iKirsch et al.l (l2005h for an 
overview. Not all of these measurements are suitable for our 
purpose. Most instruments operating only in the energy band 
below 10 keV have hard times in determining reliable spec- 
tral indices. This is because at low energies the effects of 1) 
low spectral resolution, 2) complicated instrumental effective 
area structures like filter edges etc., 3) a complex redistribution 
function and 4) the strong impact of interstellar absorption all 
conspire to make an accurate flux measurement at low energies 
and thereby a reliable photon index in the 1-10 keV band a 
hard to determine quantity. At higher energies all these effects 
are less important. For this reason, the data set that we use for 
the spectral parameters of Crab is taken from the ana l ysis o f 
BeppoSAX and other data as presented by Kui per et al.l (1200 lb . 



These are based on a careful analysis of the spectrum in a very 
broad band, separating the pulsar and nebula component, and 
although there still is a 10 % uncertainty in the absolute flux 
scale of those data, the spectral shape is very well determined, 
as elaborated below. 
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Fig. 2. BeppoSAX measurements of the Crab nebula from 



Kuiper et al.l (1200 Ih . The pulsar has been subtracted and the 
fluxes have been corrected for Galactic absorption as described 
in the text. LECS, HPGSPC and PDS data have been divided 
by 0.93, 1.01 and 0.87, respectively. For clarity, error bars (all 
lower than the plot symbols) have been omitted. The solid 
Une indicates the best fit power-law excluding LECS) with 
r = 2.147 + 0.002 and A? = 9.31 + 0.05. 



2.2. The X-ray spectrum of the Crab pulsar 



Kuiper et al. I (1200 1[ KOI hereafter) have studied Crab exten- 
sively using data from BeppoSAX, Comptel and EGRET. From 
their paper, we find the best spectrum of the pulsar to be given 
by 



Fp = 726£-'-2^^e-""^4^Vl464£-''«e-°'5^-+2021£-2.022^ 



(1) 



where E is the photon energy in keV, x = InE and Fp(E) the 
spectrum in photons m"^ s"' keV"'. This spectrum needs to be 
correcte d for absorption by the interv ening ISM, for which KOI 
used Mo rrison & McCammonI (Il983h cross sections with A^h = 
3.61 X 10^5 m-2. 



2.3. The neutron star (unpulsed) emission 

The spectrum of the pulsar as described by KOI only con- 
tains the pul s ed em ission. There is also steady emission. 
Tennant et al.l ( 1200 lb estimate the zeroth order count rate of 
the steady neutron star emission as measured by the Chandra 
LETGS as 0.19 counts s ^ For comparison, the total (dead- 
time corrected) zeroth order count rate is 8.16 counts s"\ 
hence the steady emiss ion is only 2.3 % of the total emission. 
Weisskopf et al.l (120041) estimated the steady neutron star emis- 
sion from the LETGS spectrum. For a distance of 2 kpc and 
a radius of 15.6 km, they find a 2cr upper limit to the temper- 
ature of 0.159 keV. Such a spectrum produces only 0.3-1 % 
of the total point source emission over the full LETGS range, 
1 % being reached around 0.8 keV. Hence we conclude that the 
contribution of the unpulsed emission can be neglected for the 
present purpose. 



2.4. The X-ray spectrum of the Crab nebula 

A simple power-law approximation for the nebula or the pulsar 
over the full X-ray and y-ray energy range is inadequate. At 
higher energies, the spectrum of C rab softens. This was sh own 
by KOI using Comptel data and bv lLing & WheatonI (l2003h us- 
ing BATSE data. The break of about 0.3 in photon index occurs 
around 100 keV. Ling & Wheaton give an average photon index 
of 2.1 ± 0.1 below the break at 120±40 keV and F = 2.4 + 0.2 
above the break. Therefore, we will limit the range of validity 
of any power-law approximation to £ < 100 keV. 

KOI also studied the spectrum of the nebula 
only (pulsar subtracted). These data are available at 
|http://www.sron.nl/divisions/hea/kuiper/data.html From 
their data, we find the best fit power-law based upon 
BeppoSAX MFCS data only to be F = 2.146 + 0.003 and 
- 9.31 ± 0.05 (no systematic uncertainty included). We have 
fitted all BeppoSAX data as given by KOI (MFCS, HPGSPC 
and PDS, excluding LECS) in the 1-100 keV band using a 
simple power-law, and scaling the model for HPGSPC by a 
factor of 1.01 and PDS by a factor of 0.87 with respect to 
MFCS, following KOI. We excluded LECS because it is more 
sensitive to Galactic absorption and also because its spatial 
and spectral resolution decreases rapidly for lower energies, 
making it more prone to possible systematic uncertainties. The 
best fit power-law has F = 2.147 ± 0.002 and = 9.31 ± 0.05, 
in excellent agreement with the fit to MFCS only. We show 
the fit in Fig. |2] Because the residuals are dominated by small 
remaining systematic variations of less than 1, 2 and 4 % for 
MFCS, HPGSPC and PDS, respectively, we have used an 
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Fig. 3. Comparison between the model spectra for the full Crab 
nebula and pulsar. Shown is the ratio between the power-law 
model of Toor& Seward {1914) and the more complicated 
model of iKuiper et all (l200ll) . 



unweighted fit through these data points. Hence our adopted 
model for the nebula is given by 



F„=9.31xl04£-2.i47 



(2) 



where again E is the photon energy in keV and F„(E) the neb- 
ular spectrum in photons m"^ s"' keV"' . 

In Fig. |3]we compare the model spectra of Toor & Seward 
with the combined nebula + pulsar spectrum of KOI. It is seen 
that in the MECS band (1-10 keV) both models agree within 
2 %; the differences are higher near 100 keV (up to 10 %) and 
at low energies (8 % at 0. 1 keV). 

2.5. Intrinsic curvature of the Crab nebula spectrum 

The KOI power-law model for the nebula that we use here is 
reasonably well established for energies above 2 ke V. However, 
below 2 keV interstellar absorption becomes important. The 
usual approach with low-resolution instruments has been to 
assume that the power-law extends unbroken down to about 
0.5 keV, and to adjust the interstellar column density until a 
good spectral fit at low energies is obtained. This approach was 
also used by KOI. 

There is no guarantee, however, that the power-law is un- 
broken below 2 keV. The evidence that the spectrum is not a 
simple power-law comes from both theoretical considerations 
and observations. 

2.5.1. Theoretical considerations 

The Crab nebula emits synchrotron emission due to the mo- 
tion of accelerated particles in a magnetic field. By its nature, 
such radiation has a high-energy cut-off" correspo nding to the 
maximum energy of the particles. iDel Zanna et al. have 
made detailed magneto-hydrodynamical models for pulsar neb- 
ulae with synchrotron emission incorporated self-consistently. 
Their model A was chosen to have parameters relevant for the 
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Fig. 4. Unabsorbed continuum spectrum of the Crab nebula. 
Thin solid line (F = 2.147): power-law model of iKuiper et al 



(120011) . Thick solid line: integrated spectrum derived from 
the Chandra maps o f photon index and count rate as given 
by ' Mori et al.1 (12004 ). Dashed line (F = 2.041): power-law 



approximation to the above integrated Chandra spectrum in 
the 2 -8 keV ban d . Dott ed line (F = 2.10): direct power-law 



fit of iMori et al.l ( 120041) to the full nebula spectrum. Dash- 
dotted line: theoretical model spectrum (run A) determined by 
Del Zanna et"aD(l2006h . 



Crab nebula. That model describes the spectral morphology of 
the nebula qualitatively well. The integrated spectrum shows 
multiple breaks due to spatially varying maximum particle en- 
ergies. In the X-ray band, their integrated model spectrum has 
a break around 3 keV, with a slope of 1.8 below it and 2.15 
above it (dash-dotted line in Fig.|4|l. However, this model over- 
predicts the observed X-ray flux at 1 keV by a factor of 2, and 
another model B - although less consistent with the morphol- 
ogy of the Crab nebula - has spectral breaks at different ener- 
gies. What we gather from these models is therefore that spec- 
tral breaks in the integrated X-ray spectrum may occur between 
0.5 and 10 keV, but the details of where they occur and how 
strong they are are hard to predict exactly. 

2.5.2. Observational evidence of spectral curvature 

Evidence of spectral breaks also comes from spectral imaging 
studies of the Crab nebula. There are clear spatial variations 
of the photon index over the surface of the nebula, as for ex - 
ample observed with XMM -Newton dWilUngale et al.ll200lh . 
Using Chandra ACIS-S data. lMori et all (l2004 M04 hereafter^ 
fitted simple power-laws to the local spectrum of small regions 
and produced a map of the total count rate and photon index F. 
They fixed the Galactic absorption to A^h = 3.2 x 10^^ m~^. The 
photon index ranges between 1.5 and 3.1, implying that the in- 
tegrated spectrum must show curvature. We have converted the 
0.5-8 keV count rate map of M04 into flux using version 3.7 of 
PIMMS for calculating the photon-index dependent conversion 
curve. The total nebular spectrum was calculated by adding the 
individual spectra of the small regions. This total spectrum in- 
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deed shows curvature (thick soHd Hne in Fig.|4]i, ahhough over 
the 0.3-8 keV band it never deviates more than 5 % from the 
simple power-law with F = 2.10 as determined directly from a 
spectral fit to the integrated spectrum (M04). 

Our aim is to determine how much the low-energy contin- 
uum is enhanced relative to the high-energy power-law spec- 
trum. We saw that the 1-100 keV integrated nebula spectrum 
as determined by other instruments shows no breaks (KOI; 
Sect. 12.4b . KOI put some question marks to the robustness of 
their fit below 1 keV, because of the complexity of the spec- 
tral/spatial redistribution function of the BeppoSAX LECS in 
that band. Moreover, below 2 keV the interstellar absorption 
affects the spectrum more than 2 % (below 1 keV more than 
8 %) introducing additional uncertainties. Therefore we deter- 
mine the spectral curvature relative to the well-determined 2- 
8 keV continuum. 

The spatially integrated curved spectrum in the 2-8 keV 
band is approximated by a single power-law with T = 2.041, to 
within 1 % accuracy (see Fig.|4|i. We define the intrinsic curva- 
ture correction function fdE) as the ratio of the full integrated 
spectrum to this 2-8 keV power-law approximation. The cor- 
rection factor fc{E) is shown in Fig. |5] At £ = 0.5 keV, the 
integrated spectrum is 13 % higher than the power-law approx- 
imation, showing the importance of the correction. We will use 
a simple analytical approximation for fdE): 

fc(E) = 1 + 0.41 IQ-E/0A6 (3-) 

where E is the photon energy in keV. The accuracy of this ap- 
proximation is better than 1 %, and over the RGS band (5- 
38 A) it is even better than 0.6 %. 

The curvature correction derived above (O represents the 
steepening of the soft X-rays with respect to the hard X-ray 
power-law spectrum, for the full Crab spectrum. The photon 
index of the nebula at high energies (2-8 keV) based on the 
Chandra data (2.041) differs by 0.106 from the photon index 
that we obtained in Sect. 12.41 based on KOI (2.147, see also 
Fig. nil. This difference may look surprisingly high, but can 
be explained well. The Chandra data were fitted in a low en- 
ergy band (below 8 keV), where Galactic absorption cannot be 
neglected; moreover, the Galactic column density obtained by 
M04 is 4.1 X 10^"* nr^ lower than the value used by KOI. This 
leads to different ISM transmission factors that can be approxi- 
mated to lowest order by a power-law over the 0.3-8 ke V band 
with a slope of 0.087. Thus if M04 would have used a higher 
Galactic column density (like the KOI value), their photon in- 
dices would have been higher by 0.087, close to and within the 
systematic uncertainties of the KOI value. The question then 
arises why M04 arrived at a lower column density. This may be 
due to the spectral softening of the intrinsic nebula continuum 
associated with the curvature correction; furthermore, M04 de- 
rived their column density from the fainter regions of Crab, 
which are generally in the outer parts with intrinsically steeper 
spectra. 

As the curvature correction gives the relative change in 
spectral index between the soft and hard band, and is based 
on the spatial variations of the local photon index, an overall 
change of all local photon indices due to accounting for system- 
atic uncertainties does, to first order, only affect the integrated 
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Fig. 5. Solid line: correction factor fc(E) for the intrinsic curva- 
ture of the continuum spectrum of the Crab nebula, obtained by 
comparing the full integrated model spectrum to a power-law 
approximation in the 2-8 keV band. Dashed line: analytical ap- 
proximation to fciE) (Eqn.[3). 

hard and soft band spectral slopes, but not their difference. In a 
similar way, a different treatment of the interstellar absorption 
will affect all spectral slopes in the same way and will not af- 
fect the curvature correction. Therefore, we apply the curvature 
correction (O also to the KOI continuum. 

As an independent test, we have used a MOS data set 
to derive the curvature correction in a similar way as de- 
scribed above for Chandra (we use the MOS observation of 
Appendix lA.7l ). Within the statistical errors, we obtain the same 
curvature correction. See Appendix IB. 21 for more details and 
estimates of the uncertainties on the curvature correction. 

The total, unabsorbed continuum spectrum of Crab (nebula 
plus pulsar) is then given by 

Fp(E) + F„(E)ME), (4) 

where Fp(E), F„{E) and fc{E) are given by ([T]l, (|2]l and (O, 
respectively. The systematic uncertainty on the slope of Crab 
at high energies {Fp{E) + Fn{E)) is estimated in Sect. lB.ll to be 
0.016. The systematic uncertainties on the curvature correction 
fc(E) are estimated in Sect. IB. 21 and give rise to systematic 
slope uncertainties lower than 0.010 over the RGS band. 

In order to show the effects of this curvature correction, in 
our analysis of the RGS spectrum we also show fit results for 
the case of the KOI unbroken power-law without the curvature 
correction given by ([3]l. 

3. Observations and data analysis 

3.1. Data selection 

The data were taken on February 24 and 25, 2006, during or- 
bit 1138 of XMM-Newton. We only discuss here the RGS 
observations. For each observation, either 1, 2, 4 or 8 CCDs 
were read out in a fixed sequence. This setup was done in 
order to study the effects of pile-up on the spectrum (see 
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Table 1. Observation log. 



3.2. Spectral fitting 



ID" 


seq. 


RGS 


CCDs 


^exp 


AS* 




nr 


1 or 2 


in use 


(ks) 


(arcmin) 


101 


002 




1 


15.5 


- 


101 


003 




2 


15.2 


- 


101 


004 




3 


7.0 


- 


101 


005 




4 


3.7 


- 


101 


006 




5 


1.6 


- 


101 


007 




6 


0.6 


- 


101 


008 




8 


0.5 


- 


101 


009 




9 


1.5 


- 


101 


033 


2 


1 


15.5 


- 


101 


034 


2 


2 


15.2 


- 


101 


035 


2 


3 


7.0 


- 


101 


036 


2 


5 


3.7 


- 


101 


037 


2 


6 


1.6 


- 


101 


038 


2 


7 


0.6 


- 


101 


039 


2 


8 


0.5 


- 


101 


040 


2 


9 


1.5 


- 


101 


010 


1 


1,5 


4.0 


- 


101 


Oil 


1 


2,6 


4.5 


- 


101 


012 


1 


3,8 


4.0 


- 


101 


013 


1 


4,9 


4.5 


- 


101 


041 


2 


1,6 


4.0 


- 


101 


042 


2 


2,7 


4.5 


- 


101 


043 


2 


3,8 


4.0 


- 


101 


044 


2 


5,9 


4.5 


- 


101 


014 


1 


1,3,5,9 


3.9 


- 


101 


015 


1 


2,4,6,8 


2.5 


- 


101 


045 


2 


1,3,5,7 


2.5 




101 


046 


2 


2,6,8,9 


2.9 




101 


016 


1 


all 


1.7 




101 


047 


2 


all 


1.9 




201 


004 


1 


all 


5.8 


+3.5 


201 


005 


2 


all 


5.8 


+3.5 


301 


004 


1 


all 


5.7 


-3.5 


301 


005 


2 


all 


5.7 


-3.5 


401 


004 


1 


all 


8.5 


+7.0 


401 


005 


2 


all 


8.5 


+7.0 


501 


004 


1 


all 


14.9 


-7.0 


501 


005 


2 


all 


14.9 


-7.0 



" Observation ID (the leading digits 0312790 are the same for all ob- 
servations and have been omitted here) 
* off-axis angle in the cross-dispersion direction 



Appendix lA.4b . In addition, a set of observations with off-axis 
pointing in the cross-dispersion direction was done, with off- 
axis angles of ±3.5 and ±7.0 arcmin, in order to study the 
spatial distribution of the emission in the cross-dispersion di- 
rection (see Appendix IA.6b . A summary of the observations 
is shown in Table [T] The observations were processed using 
SAS version 7.0 and the set of calibration files (CCF) avail- 
able October 2006. More information on the spectral extrac- 
tion and the various correction factors that we need are given 
in Appendix A. 



All spectral fitting was done using the SPEX package 
dKaastra et al.l Il99a, see also www.sron.nl/spex). RGSl and 
RGS2 data were fit simultaneously. The continuum model for 
Crab was taken to be the KOI model (pulsar plus nebula) dis- 
cussed in Sect. 12.21 and 12.41 including the curvature correction 
derived in Sect. 12.5.21 This continuum model was read as a 
table (model file in SPEX), and the normalisation of this con- 
tinuum was frozen. We apply two multiplicative components to 
this continuum, discussed below: a correction for the unknown 
effective area of the RGS, and a model for the interstellar ab- 
sorption. Therefore in our fit, the only free parameters are the 
properties of the interstellar absorption (column densities and 
line width), and the two parameters of the effective area correc- 
tion. 

As we indicated in Appendix |A.2| the effective area of RGS 
has been determined modulo a power-law correction. The true 
effective area At„s(A) as a function of the wavelength A (in A) 
can be related to the effective area currently present in the re- 
sponse matrix Aresp(/l) by AtrueC/l) - fiA)Aresp{A) where 



/(^)=A(^/10)". 



(5) 



For technical reasons we have chosen not to fit directly the pa- 
rameters A and a, but to fit use f(Ai) and /(/I2) with /1 12 ex- 
treme points for the RGS band, chosen here as 5 and 40 A. 

For our absorption model, we use a combination of the hot 
and slab models of SPEX. Both models use the same set of 
atomic absorption cross sections (see below). The slab model 
calculates the transmission of a layer of plasma with arbitrary 
composition. Free parameters are the column densities of indi- 
vidual ions as well as the intrinsic velocity broadening. We use 
this model component for all neutral atoms of the elements that 
show clearly detectable absorption edges in the RGS band (N, 
O, Ne, Mg and Fe), plus their singly ionised ions (for detections 
or upper limits to the ionised fraction). 

For all other elements we use the hot model. This model 
calculates the transmission of a plasma in collisional ionisation 
equilibrium. We use this with a very low temperature (the neu- 
tral gas limit). For the abunda nces of those el ements we adopt 
the proto-solar abundances of iLoddersI (120031) . The dominant 
contributors in the RGS band for those elements are He, H, C, 
Si and S (in order of importance). The cross sections of all these 
elements are tied to the hydrogen cross section. 

The cross sections used for both the hot and slab models 
are described in more detail in the SPEX manual. Basically, 
both continuum and line opacity are taken into account. For 
neutral atoms in the RGS band predominantly inner shell ab- 
sorption lines are important; the cross sections of those line 
were mainly calculate d by E. Behar using t he HULLAC code; 
see also Section 5.1 of Kaastra et al.l (120021) for more details. 

Alternatively, in order to investigate the role of dust, we 
have also made spectral fits with the model outlined above 
(which we label model A), but with a part of the total col- 
umn replaced by a gas and dust mixture (model B). This gas 
and dust mixture is assumed to have also the proto-solar abun- 
dances of LoddersI (12003) . but a fraction of the gas is depleted 
into dust grains. The depletion factors of the dust are taken from 
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Wilms et al. I (l2000l) . The transmission of the dust is calculated 
using the dabs model that we introduce here in the SPEX pack- 
ag e, which fo llows completely the dust treatment as described 
by IWilms et a l. (2000) and references therein. Basically, the 



Table 2. Best fit parameters. 



grains are assumed to be spherical, fluffy (p - 1000 kgm" ), 
and have a size distribution dn/da ~ a^'' for the grain radius a 
between < a < flmax, with default values ami,, = 0.025 fj.m, 
dmax - 0.25 //m, and p - 3.5. 

As a third model (Model C) we also made a fit without the 
curvature correction but further identical to model A, in order 
to investigate the effects of this curvature correction. 

For the spectral fitting, we rebinned all data by a factor of 
5 to account for the reduced spectral resolution caused by the 
spatial extent of Crab (0.22 A FWHM, see AppendixlATli. This 
gives bins with a width of ~ 0.04 A at 7 A to 0.07 A at 30 A. 
Fitting was restricted in most cases to the 7-30 A band; be- 
low 7 A, scattering by the RGS gratings becomes very impor- 
tant, but theory of scattering at these short wavelengths starts 
to break down. Therefore, there is essentially no knowledge of 
the RGS effective area below 7 A. Hence this wavelength range 
was excluded from the fit. Above 30 A, the low energy of the 
X-ray events deposited on the CCD, combined with the extent 
of the source, starts to coincide with the tail of the CCD system 
peak. For this reason, PI selection was limited at low energies, 
as explained in Appendix lA.il Still, some contamination of X- 
ray events by the CCD system peak can be expected. For this 
reason, the spectrum above 30 A was also excluded from the 
spectral fit (see Fig. |6]l. However, we were able to derive the 
nitrogen abundance by extending the fit range - for the sole 
purpose of determining that abundance - to 32 A. That fit gives 
biased values for all other parameters, but the jump across the 
nitrogen edge at 30.8 A is determined well. 

Furthermore, we omitted small regions near the interstel- 
lar oxygen K-edge (22.7-23.2 A) and iron L-edge (17.2- 
17.7 A), for reasons explained in more detail in our discussion 
(Sect. 14. 1 .2|j4T3l essentially the presence of molecules that 
affect the fine structure of the edge). As the ISM abundances 
are mainly determined from the depth of the absorption edges, 
ignoring these small regions does not affect the basic outcome 
of our analysis. 



parameter" 


model A 


model B 


model C 




pure gas 


gas & dust 


pure gas 




curv. corr. 


curv. corr 


no curv. corr 


X' 


2341 


2327 


2366 


d.o.f. 


1316 


1315 


1316 




1.056 ±0.009 


1.055 ±0.009 


1.035 ±0.009 


/(40)'' 


0.746 ± 0.016 


0.762 ±0.016 


0.854 ±0.019 


column densities'': 


Hi'''' 


31750 ±220 


9700 ± 1700 


31280 ±220 


Hi''./ 




22300 ± 1800 






31750 ±220 


32060 ± 220 


31280 ±220 


Hii'' 


1752 ±4 


1752 ±4 


1752 ±4 


Ni 


2.67 ± 0.24 


2.72 ± 0.24 


2.48 ± 0.25 


Oi 


19.28 ±0.18 


19.83 ±0.22 


19.23 ±0.18 


On 


0.23 ± 0.14 


0.25 ±0.15 


0.20 ±0.14 


Nei 


4.69 ± 0.21 


4.66 ±0.21 


4.74 ± 0.21 


Neil 


0.46 ± 0.25 


0.52 ± 0.24 


0.64 ± 0.25 


Mgi 


1.14 ±0.28 


1.20 ± 0.28 


1.13 ±0.28 


Mgn 


0.00 ± 0.09 


0.00 ± 0.10 


0.00 ± 0.09 


Fei 


0.65 ± 0.05 


0.77 ± 0.04 


0.65 ± 0.05 


Fell 


0.23 ± 0.06 


0.14 ±0.04 


0.26 ± 0.07 



" Error bars are all Icr for one interesting parameter. 

* /(5) and /(40) are the values of the effective area scale factor /(/I), 

Eqn. (O, at 5 and 40 A. 

'■ Ionic column densities are given in units of 10^' vcT^. The column 
densities of all metals include all phases (dust and gas). 
'' H I includes in the cross section contributions from neutral He, C, Si, 
S, Ar, Ca and Ni assuming proto-solar abundances for these elements. 

For the pure gaseous phase only. 
/ For the gas & dust mixture with proto-solar abundances added to 
the pure gaseous phase above. 

^ Total column: sum of contributions listed as b and c. 

'' For Hn the column is derived from the pulsar dispersion measure 

assuming that all free electrons originate from hydrogen (see text). 



index p. However, our fit is not very sensitive to these param- 
eters, and therefore we have kept them to their default values 
listed in Sect.[l2l 



3.3. Results of the spectral fits 

We have obtained a fit using model A outlined in Sect. 13.21 
with ax^ of 2341 for 1316 degrees of freedom. The fit is shown 
in Fig.|6]and the parameters are given in Table|2l The fit residu- 
als over the fitted range have - on average - an intrinsic system- 
atic scatter of 1.8 % and a random noise component of 1.9 %, 
yielding a total scatter of 2.6 %. Therefore, the typical system- 
atic uncertainty in our model is less than 2 %. The higher devi- 
ations outside the fitted range can be simply understood as we 
already discussed in Sect. 13.21 

For our alternative model B that includes dust, we list the 
best fit parameters also in Table|2] We have tried to leave the ba- 
sic dust parameters free in our fits: grain mass density p, mini- 
mum and maximum size a^oin and Cmax, and the size distribution 



3.3.1 . Search for hot gas 

The fit residuals near 21.6 A (Fig.|7]l and 19.0 A (Fig.|7]l, where 
we would expect the 1 s-2p transitions of O vii and O viii re- 
spectively, show some evidence of additional weak absorption 
features. We have determined the formal equivalent width of 
these features by adding narrow absorption lines to the model 
(but broadened due to the spatial extent of Crab as seen by 
RGS). 

For the O vii line, we find a formal equivalent width of 9 ± 
4 mA, and for O viii 7 + 2 mA if the line is narrow. For O viii the 
feature looks somewhat broader; if we allow the width of the 
line to vary, we find a higher equivalent width of 18 ±4 mA, but 
this requires a high velocity width (Gaussian cr - 5000 km s '). 

For consistency, we have checked a Chandra LETGS spec- 
trum of the Crab pulsar (Obsid 759), and found a formal equiv- 
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20 

Wavelength (A) 

Fig. 6. Best fit RGS spectrum of the Crab nebula using model A. The fit was restricted to the 7-30 A range (indicated by the 
vertical lines), excluding two small regions near the O i K-edge and Fe i L-edge (also indicated by vertical lines). However in the 
plot we show the full data range. Single and double CCD readout modes have been combined, but RGSl (dark colour) and RGS2 
(light colour) data are shown separately. The upper panel shows the fluxed RGS spectrum (but at the original spectral resolution 
of the data, including spectral broadening due to the spatial extent of the nebula). For clarity reasons we plot the flux in keV units 
along the y-axis, which reduces the dynamical range along that axis significantly. The lower panel shows the relative fit residuals. 



alent width of the O vii line of 7 ± 17 mA, consistent with our 4.1 .1 . Neon K-edge 
RGS measurement. 



4. Discussion of the absorption spectrum of Crab 

4.1. Fine structure near absorption edges 

Our fits (Fig.|6]) clearly show strong absorption edges of neutral 
O, Fe and Ne. Other edges are weaker (like Mg at 9.47 A and N 
at 30.77 A). In principle, the neutral Si edge at 6.72 A is also in 
the RGS band and is expected to have an optical depth of 0.017; 
however, it is in the part of the spectrum where our fits be- 
come poor and also any uncertainty related to the instrumental 
Si edge will contribute to its systematic uncertainty. Therefore 
we will not discuss Si further in our paper. We will discuss here 
in more detail the edges of Ne, Fe, and O. 



The Ne-edge is clearly detected in the data (Fig. |7]i. The opti- 
cal depth T of the edge is 0.15. Nearby is the Li edge of neu- 
tral iron, but that edge is much weaker: r - 0.010. Recently, 



Juett et al. ( 2006 ) have studied the fine structure near the neon 



K-edge and iron L-edge using Chandra HETGS spectra of 
X-ray binaries. The spectral resolution of the HETGS (0.01- 
0.02 A), combined with the point source nature of these ob- 
jects, gives them effectively a 10-20 times higher spectral res- 
olution than our Crab spectrum (extended source). For Nei, 
in addition to the pure absorption edge, there is the series of 
Is-np transitions (n > 2 because the 2p shell is completely 
filled for Ne i). According to the model of Juett et al., the ls-3p 
line is at 14.295 + 0.003 A, and this line is strong enough that 
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17 IE 
Wavelength (A) 




22 23 
Wavelength (A) 

Fig. 7. Blow-up of Fig. |6]near the Ne i K-edge (top), Fe i 
edge (middle) and Oi K-edge (bottom). Crosses: RGSl; dia- 
monds: RGS2; solid line: model A (no dust); dotted line: model 
B (with dust). 



the effective K-edge of neon shifts towards this wavelength. 



For the true edg^j, we used the value of 14.249 A. When we 
leave it as a free parameter, the best fit value for the edge is 
14.231 + 0.013 A. This low shift is most likely caused by our 
neglect of the Rydberg series for n > 5 in our absorption model. 

The K-edges of O viii and Ne i differ only by ~ 0.02 A. With 
our current effective spectral resolution for Crab, we cannot 
resolve these edges. However, as we show later (Sect. 14.3. Il l 
the amount of hot gas in the line of sight is small enough to be 
neglected. 

4.1 .2. Iron L-edge 

We have already seen that the Li edge of Fe i is very weak, and 
we will not discuss that edge further The predicted depth of 
the L2 and L3 edges for our best-fit model A are only t - 0.027 
and T - 0.056, respectively. At a first glance this looks con- 
tradictory to the deep features seen in Fig.|2] The difference is 
due to the presence of strong inner-shell line transitions (for 
example 2p-3d excitations). Such tra nsitions have first been 
identified in an astrophysical source bv lSako et al.l (1200 ll) . and 
have subsequently been incorporated in our absorption mod- 
els in SPEX. Those transitions cause the two deep absorption 

troughs at 17.1 and 17.4 A. They are also visible in the HETGS 

I 1 J — — I 

data of lJuett et al.l (12006), at similar wavelengths. At our ef- 
fective spectral resolution it is not possible to separate the in- 
ner shell transitions of Fe i from those of Fe n; the blends for 
these ions almost coincide. Therefore, the difference between 
the absorption by both ions is mainly derived from the different 
location of the absorption edges (see Fig.|7]i. 

Similar to the case of the neon edge, also here possible 
contamination by highly ionised oxygen must be considered. 
Within our effective resolution, the O vn K-edge is close to the 
Fe II L3 edge (Fig.|7]i. Our model gives an optical depth of 0.019 
for this Fe 11 L3 edge, while the O vii edge has a depth of at most 
0.002. Thus, contamination of the iron edge by highly ionised 
oxygen is less than 10 %, well below the statistical uncertainty 
on the Fe 11 column density. 

Although the column density of Fe is determined well from 
the jump across the absorption edge, our fit is not optimal. At 
the long wavelength side of the edge, between 17.5-18.0 A, 
there is excess flux of the order of 5 %. This excess extends 
downwards to 17.1 A. Adopting a ~ 5 % higher continuum 
between 17-18 A, our absorption model provides a reasonable 
fit. Note that our model B with dust provides a slightly better fit 
in the absorption lines compared to model A; this is due to the 
high opacity of the lines, which results into stronger shielding 
effects by dust grains, and therefore to a relatively lower line to 
continuum opacity than for the pure atomic case. 

The 5 % amplitude of the excess is much higher than any 
uncertainties in our cross-dispersion correction factor, exclud- 
ing therefore an origin in uncertainties of dust scattering near 
the iron edge. As Crab has a rather smooth synchrotron contin- 
uum, we do not expect such broad features to be present in the 
intrinsic Crab spectrum. 



' The definition of "edges" by Juett et al. is non-standard. They de- 
fine the edge where empirically the flux drops; we keep here the more 
formal definition of the series limit of the absorption lines. 
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It i s striking, however, that several of the fits of iJuett et al 
(l2006h to spectra of X-ray binaries have excess flux of a similar 
magnitude just at the long wavelength side of the Fe edge, sim- 
ilar to our case (for example Cyg X-1). Although their features 
look a little sharper, convolving them with the spatial/spectral 
response of our Crab observation would yield similar features. 

The instrumental fluorine K-edge near 17.8 A is determined 
well from higher resolution RGS spectra of point sources with 
low interstellar absorption columns. Those spectra do not show 
similar residuals. Hence we can exclude an origin of the excess 
in the instrumental F K-edge. 

Finally we note that it is unlikely t hat this excess is du e 
to line emission from a supernova shell. ISeward et ak ( 2006h . 
using Chandra data, have derived an upper hmit to the lumi- 
nosity of such a shell of 2 x 10^^ W. The observed excess 
with RGS, when fitted with a Gaussian profile, has a flux of 
109 + 8 phm"^ s at a mean wavelength of 17.53 + 0.03 A and 
with cr = 0.57 ± 0.08 A. At the distance of Crab, and correcting 
for foreground absorption, this corresponds to a (line) luminos- 
ity of 2.4 X 10^^ W, almost equal to the upper limit to the total 
X-ray luminosity of any shell. Moreover, RGS cannot see all 
flux from such a shell unless the radius is less than 2.5'. For 
such a radius, the maximum wavelength shift observed with 
RGS due to the spatial extent of the emission is 0.35 A. There 
are no strong line transitions close enough to explain the emis- 
sion. For example, if the emission would be O viii Lya emission 
at 18.97 A, it should be located about 10' from the Crab pulsar. 

4.1 .3. Oxygen K-edge 

The oxygen K-edge clearly is the strongest feature of the spec- 
trum, with a jump of a factor of 2 across the edge. The Is- 
2p absorption line is reproduced well by our fit. However, in 
our original fit the edge energy, for which we used a value of 
23.045 A, clearly seems to be off. Making it a free parameter, 
we find a best-fit wavelength of 22.917 ± 0.008 A, a shift of 
-0.128 A. This does not mean that the edge itself is shifted. 
Two effects may cause the discrepancy. 

First, we use the continuum cross sections of 



Verner & Yakovlevi (Il995h . In this calculation, there is no fine 
structure at the edge, which is located at 538 eV (23.045 A). 
The edge, however, has at least two strong sub-edges at the 
2f series limit (22.567 A), and series limit (22.767 A). 
Furthermore, the calculated wavelengths of the absorption 
lines that we use cannot be calculated more accurately than 
typically 10-30 mA. Only the strongest transition (ls-2p) 
was adjusted empirically to observed values. All this causes 
significant uncertainty in the adopted absorption model for 
atomic oxygen. We do not attempt to resolve th is problem here 
but refer to a forthcoming paper dde Vries et al. 2008) for a 
more extensive discussion of the K-shell absorp tion structure 
in X-r ay sources; al so see the discussion in iPaerels et al 
(1200 lb and lJuett et all ([2004). 



Secondly, we note that in our fit with model B, with our 
adopted oxygen dust depletion factor of 0.4 for the dust+gas 
mixture, the total fraction of oxygen bound in dust is 28 % 
(as there is also a pure gaseous component). This number is 



quite uncertain as it depends on our adopted dust depletion 
factor Oxygen bound in minerals or dust grains has a clearly 
different fine structure near the e dge than atomic oxyge n; see 
for instance the many examples in Ivan Aken et al.l ( Il998h . This 
changes the absorption structure for the chemically bound oxy- 
gen population significantly. 

As the O II and O iii edges are separated well from the O i 
edge, the high spectral resolution of RGS allows us to put tight 
upper limits on the amount of ionised oxygen. No strong edges 
of these ions are seen in the spectrum. 

Fig.|2]also shows a small depression near 21.6 A, the wave- 
length of the strong resonance line of O vii. We discuss that 
furflier in Sect. 14. 3. II 

Fin ally, based on our wo rk on blazars as described in pa- 
per III dKaastra et al. l2008ah . we estimate that the systematic 
uncertainty on the neutral oxygen column in Crab due to the 
uncertainty in the instrumental oxygen edge is less than 1 %. 

4.2. Column densities of tine cold gas 

From the column densities obtained from Table |2] we see that 
the absorbing medium is mostly neutral, except for iron, for 
which we find a significant ionised contribution. Column den- 
sities towards Crab ha ve also been measured using UV and 
optical spectroscopy by iSollerman et al] (l200d SOO hereafter). 
They derived column densities based on a Doppler width of 
1 kms For such a low velocity broadening, the lines are 
strongly saturated and hence the statistical uncertainty on the 
column densities are relatively high, despite significant line de- 
tections. 

SOO obtained a neutral hydrogen column density of (3.0 ± 
0.5) X 10^^ m"^, fully consistent with our more precise value of 
(3 . 1 8 ± 0.02) X 1 0^5 m-2 for model A or (3 .2 1 + 0.02) x 1 0^^ m'^ 
for model B. 

There are 4 ions in common between our work and SOO: 
O I, Mg I, Mg II and Fe ii. For O i and Fe ii, the column densities 
agree within the statistical errors of a factor of 2 and 4, respec- 
tively. However, for Mg i and Mg ii there is a bigger difference. 
SOO find a Mg ii/Mg i ratio of 5, with an uncertainty of a factor 
of 2, while we find a ratio less than 0. 1 . The total magnesium 
columns however agree: we find (1.1 - 1.2 + 0.3) x 10^^ m"^ for 
models A and B respectively, compared to (0.7^^4) x 10^^ m"^ 
for SOO. 

The different ratio appears puzzling but can be understood 
as follows. First, the Mg i and Mg 11 edges are close, at 9.48 and 
9.30 A. The optical depth of the Mg i edge (our deepest Mg 
edge) is about 0.02. In this region of the spectrum (see Fig.|6]l, 
there is a small difference between the RGS 1 and RGS2 data, at 
the 1 % level. It is therefore easy to get systematic errors at this 
level. While the total magnesium column is rather robust (as it 
measures essentially the jump across the combined nearby Mg i 
and Mg 11 edges), discriminating between both ions is harder 
because this depends solely on the small data stretch between 
9.30-9.48 A. 

We therefore believe that a higher Mg ii/Mg i ratio than re- 
ported by us is quite well possible, but that our total magnesium 
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column is rather secure. Such a higher ratio would also be con- 
sistent with the enhanced degree of ionisation of iron. 

4.3. ISM constituents 

For a proper discussion of abundances, the different phases of 
the ISM should be taken into account. Near the Sun, we have 
th e following com ponents (see the recent overview of the ISM 



f K 

bv lFerrierel ( 120011) for more details): 



hot ionised gas (~ 10^ K), 
warm ionised gas (~ 8000 K), 
warm atomic gas (6000-10000 K), 
cold atomic gas (20-50 K), 
molecular gas (10-20 K), 
and dust 



We consider these components in more detail below. 
4.3.1 . Hot ionised gas 

In general, the presence of hot gas in the line of sight towards 
any X-ray source is most easily seen using high spectral reso- 
lution measurements of the associated absorption lines. This is 
because in general the column densities of highly ionised gas 
are relatively low and therefore the absorption edges are weak, 
while the absorption lines still may be visible. 

X-ray absorption due to the p resence of hot gas in our 
Galaxy has been reported before. Yao & Wang (|2006) anal- 
ysed Chandra grating spectra toward the bright low mass X- 
ray binary 4U 1820-303, which has a neutral hydrogen column 
density of 1.5 x 10^^ m"^, half of the value towards Crab. In 
4U 1820-303, these authors find equivalent widths for the O vii 
and Ovm ls-2p absorption lines of 40 ± 12 and 19 + 6 mA, 
respectively. These values are much higher than the values that 
we obtain for Crab (Sect. [33Tt . of 9 + 4 and 7 + 2 mA, respec- 
tively. 

Due to the weakness of the O vii and O viii absorption lines 
toward Crab, it is hard to estimate the associated column den- 
sities. If the velocity dispersion of the gas o-y = 180 kms"', 
the value found by Yao & Wang toward 4U 1820-303, then 
the column densities of O vii and O vm are (8 + 3) X lO'** and 
(2.6 ± 0.9) X 10'^ m"^, respectively, 2.5 and 10 times below the 
values towards 4U 1820-303. The optical depth at the K-edge 
of these ions is then less than 0.002, and therefore the neglect 
of these ions in our analysis is justified. For solar abundances, 
and assuming that all oxygen of the hot gas would be in the 
ionisation stage of O vii or O vm, the total hydrogen column 
density of the hot gas would be 2 x 10^^ m"^, less than 1 % of 
the neutral hydrogen column density. 

Unfortunately, we do not have a good constraint on o-y. 
If (Ty were substantially lower than 180 kms much higher 
columns would be required. We have worked this out as fol- 
lows. We made a fit to our spectrum adding columns of Ovii, 
O VIII and Ne ix to our model in the form of an additional slab 
component. These ions are dominant X-ray opacity sources for 
typical temperatures of a million K. Surprisingly, our fit im- 
proves significantly {Ax^ = 88). The improvement is produced 



mainly by the inclusion of the absorption lines as well as fine 
tuning of the continuum, in particular near the neon and iron 
edges, although the ~ 5 % excess emission near the Fe-L edge 
still remains. Because the O vii K-edge is close to the Fe ii L2.3 
edges (see Fig.|7]l, the fit reduces the Fe 11 column to zero in ex- 
change to the O VII column. In a similar way, Ne i is reduced in 
favour of Oviii, also because both edges are so close (Fig.|7]i. 
Still the overabundance of neon remains. 

However, this model is not acceptable because of other rea- 
sons. First, the total column of highly ionised oxygen (O vii and 
Ovm) is high: 5 x 10^' m"^. Combined with the ionised hy- 
drogen column deduced from the electron density (determined 
from the pulsar dispersion measure, see Sect. 14.3.2b . an over- 
abundance of oxygen of 6 times solar is required. This over- 
abundance is even a lower limit if a significant fraction of the 
ionised hydrogen is associated with warm (10"* K like) gas in- 
stead of hot (10^ K like) gas. The most important problem with 
this model is, however, that the required velocity dispersion of 
the gas is cTy < 3 kms"', far below the thermal velocity of the 
hot gas. If the hot gas is predominantly in collisional ionisa- 
tion equilibrium, then for a typical temperature of lO*" K, the 
thermal velocity of oxygen ions is more than 20 km s"' . 

Of course, we can solve this problem by accounting prop- 
erly for the thermal width of the lines. We have done this by 
replacing the additional slab component by a hot component 
with solar abundances. Additional free parameters are the to- 
tal hydrogen column density, the temperature and the velocity 
broadening o-y to be added in quadrature to the thermal width 
of the lines. This model yields only a modest fit improvement 
over the original model without hot gas: A^^ = 23. We find a 
typical temperature of 2-3 million K and a hydrogen column 
of 5 X lO^^* m"^, almost two orders of magnitude below the 
column density of neutral hydrogen. Accordingly, the optical 
depths of the strongest oxygen edges for this component are 
at most 0.001, and there is no significant change to any of the 
parameters for the cold absorption component. 

In summary, we conclude that some hot gas may be present 
in the line of sight towards Crab, at a level of a percent of the 
neutral hydrogen column. Typically, the amount of hot gas to- 
wards Crab is less than A^h ^ 5x10^^ m"^. This is small enough 
that neglecting it does not affect our results for the cold phase 
or the calibration of the Crab spectrum. 

As the abundances in the hot ionised gas may be slig htly 
diff'erent from those of the colder components ( Yao & W angI 
l2006l) . we exclude this hot component explicitly from our abun- 
dance analysis. 



4.3.2. Warm ionised gas 



Pulsa r dispersion measurements (for example iLundgren et al 



I995I) yield a total electron column density of ( 1 .752 + 0.003) x 



10^^ m"^, where the error is dominated by the monthly fluctua- 
tions of this number, proba bly caused by the Crab nebula itself 
( Isaacman & Rankinlll977h . As this electron column density is 
high compared to the hot ionised gas column density, most of 
the free electrons must be produced by the warm ionised gas. 
For the typical conditions prevailing in that medium, most of 
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these electrons come from hydrogen. Hence, we use the elec- 
tron column density as an estimate for the H n column density. 

The other elements in the warm ionised gas are either neu- 
tral or singly ionised, and they are accounted for in our absorp- 
tion model, although we cannot distinguish between the neutral 
atoms in this gas phase and those in the cold and warm atomic 
gas. 

4.3.3. Cold and warm atomic gas 

These two gas phases mainly contain neutral atoms. In the X- 
ray band, they are hard to distinguish, the main difference be- 
ing the thermal line width. These are hard to measure, how- 
ever, as the associated absorption lines are strongly saturated. 
Therefore we discuss the combined abundance of both phases, 
and also include the neutral and singly ionised atoms from the 
warm ionised gas in this component. 

4.3.4. Molecular gas 

Special attention needs to be paid to molecular hydrogen (H2) 
here, as it is the main molecular X-ray opacity source. 

The distribution of molecular gas in the Galaxy is very ir- 
regular and depends strongly upon the line of sight. It has a 
relatively high concentration closer to the centre of the Galaxy 
and is less important further outwards. In the local ISM, about 
20-25 % of the hydrogen is in its molecular form dSavage et al.l 
1977h . and it is this value that is used as default in for ex- 
ample the absorption model of Wil ms et al. I (I2OOO), although 
these authors warn the user to take care of such spatial varia- 
tions. Over the X-ray band, the H2 opacity per hydrogen atom is 
about cTmoi / CTatom - 1 -42 timcs higher than for H i (Wilm Tet al.l 
2OOOI) . Hence, for a mixture with a 20 % molecular contribu- 
tion, the effective hydrogen opacity should be 8 % higher than 
for pure H i; ignoring this molecular component would there- 
fore bias abundance determinations. 

In our spectral model we have ignored molecules, so we 
need to justify this. We do this by determining R, the rela- 
tive enhancement of the total opacity including molecules com- 
pared to the opacity of the same gas if it is fully atomic. R can 
be written as 



O"mol 
C^atom 



l)x X/h 



(6) 



where here n^oi - 2nu2 and «atom are the number densities of 
hydrogen nuclei in molecular and atomic form, respectively. 
Furthermore, /h is defined here as the effective opacity of hy- 
drogen relative to the other metals in our X-ray spectrum. 

For gas with solar composition, the contribution of hydro- 
gen to the total opacity varies between 5%atlOAtol7% 
at 30 A. A typical weighted average value for our Crab spec- 
trum is ~ 10 %. As approximately half of the total opacity is 
produced by elements such as N, O, Ne, Mg and Fe for which 
we determine the column densities separately, we estimate /h 
by comparing the pure hydrogen opacity to the total opacity 
produced by the elements that we tied to hydrogen using solar 
abundances. This leads to an estimated of /h - 0.2. 



Crab lies in the direction of the Galactic anticentre, at a 
distance of 2 kpc. Therefore the fraction of H2 is expected to 
be lower than t he typ ical value of 20 %. From the CO maps 
of lOame et alj ( 1200 lb (available through NASA's skyview fa- 
cility), and using their conversion factor from CO equivalent 
width to H2 column, we find nm < 5 x 10^^ m"^. There is 
a molecular cloud with a higher column density of 10^"* m"^ 
North of Crab, but it is more than 5' away. Therefore, our 



best estimate is nmoi/«at 



0.003 (direct upper limit) or 



nmoi/natom < 0.06 (worst case scenario, if the cloud North of 
Crab would cover it). 

Taking all these factors into account, we find R < 0.42 x 
0.003 X 0.2 < 3 X 10""*. For the worst case scenario with 



Wmol/na 



< 0.06, still R < 0.005. Therefore, we can com- 



pletely ignore the role of molecular gas in our line of sight. 
4.3.5. Dust 

Dust in the line of sight towards Crab or any other X-ray source 
has two effects: scattering and absorption of X-rays. 

Scattering does not destroy photons but merely redis- 
tr ibutes th em into a halo around the source. Recently, 



Seward et al. (2006) have studied the dust scattering halo of 
Crab using Chandra imaging. They conclude that the total col- 
umn density of scattering material (dust and associated gas) is 
1.7 - 2.2 X 10^^ m"2 for their best fit model, and 2/3 of that 
value for an alternative model. Most of this can be modelled to 
be in a smoothly distributed component between us and Crab, 
while 10-25 % of the scattering material is in a single struc- 
ture at a distance of ~ 100 pc from the Sun, perhaps associated 
with the boundary of the local hot bubble. Surprisingly, the to- 
tal column density that Seward et al. (2006) find is less than 
the column derived from X-ray absorption measurements. This 
may indicate a relative low abundance of dust towards Crab. 

The amount of dust inside Crab is relatively small, and can 
contribute only 1-2 % to the total column density t owards Crab. 



This can be derived from infrared observations: Temim et al 



( 2006) estimate that Crab contains about 0.004-0.010 Mq of 
graphite or 0.006-0.015 M© of silicates, in order to explain the 
thermal dust emission seen at 70 fim. Assuming a constant den- 
sity sphere with radius 2', the total carbon column of the dust 
within Crab (in case of graphite) would be 0.7 - 1 .8 x 10^" m"^, 
and in case of silicates (assuming for simplicity olivine or 
Fe2Si04) the silicon column would be 0.6 - 1.6 x 10'^ m"^. 
In both cases, this corresponds to 1-2 % of the total interstellar 
column. Hence we will ignore this component. 

Apart from scattering, dust also absorbs X-rays. There 
are two main differences between the absorption properties of 
atoms contained in dust and free atoms. First, the absorption 
edges and lines in bound atoms are slightly different, due to the 
chemical shifts. Secondly, the effective absorption cross section 
can be less due to shielding in grains. 

Unfortunately, our spectral resolution is not very high due 
to the extended nature of the Crab emission. Therefore it is 
hard to measure accurately the location of absorption lines. 
However, our data still contain some usefull information on 
that. 
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Constituent 


oxygen 


iron 


iron 


ref. 




main line 


main line 


2nd line 




Ferrous, Fe'^ 




17.498 


17.196 


a 


Olivine 


23.26 






b 


Ferric, Fe'^ 




17.456 


17.130 


a 


Hematite 


23.43 






b 


Atomic O I 


23.508 






c 


Atomic Fe i 




17.453 


17.142 


d 


Crab 


23.466 ± 0.009 


17.396 ± 0.009 


17.120 ±0.016 


e 



" Ivan Aken & Liebsched ( l2002h . corrected by -0.019 A following 
ljuett et al.l (12006 



van Aken etai](ll998h 



Juett et al. ( 2004), from Chandra spectra 
Behar et al. C2001,) , HULLAC calculations 



■ Present work 



In Table [3] we compile our current knowledge about the 
dominant absorption lines in different constituents, and com- 
pare it with the observed centroids as measured in our Crab 
spectrum. These observed wavelengths were obtained by ignor- 
ing the absorption lines in our best fit model and fitting them 
instead with narrow absorption lines. Olivine (Mg2-AFej:Si04) 
is one of the best known examples of a ferrous, and hematite 
(Fe203) of a ferric compound. 

If we consider neutral oxygen, then the observed absorp- 
tion line centroid suggests approximately equal contributions 
from atomic and ferric iron (Table|3]l. Such a mixture also gives 
good agreement with the secondary iron line blend at 17.12 A, 
and is also in agreement with the ~ 30 % of oxygen bound 
in dust derived from our model B. There is only disagreement 
with the main iron line blend at 17.40 A. Clearly, both ferric 
and atomic iron have their predicted line closer to the observed 
centroid of this blend than feiTous iron, but none of them agree 



fully. I t may be that the predicted wavelengths bv Behar et al 



( 1200 lb are slightly off, as these lack detailed experimental vah- 
dation; however, their calculations agree within a few mA with 
other calculations (Raassen & Kaastra, private communication, 
as quoted by Behar et al.). We suggest that the dust line of sight 
towards Crab is dominated by ferric compounds. 

The second difference in absorption properties of dusty 
atoms compared to free atoms is the shielding of atoms in 
dust grains. The shielding leads to effectively lower opacity per 
atom. This is illustrated in Fig. [8] At the oxygen edge, the dif- 
ference can be up to 14 %. In practice, the effects of inclusion 
of dust in the model are not as dramatic as shown here. The 
reason is that one can always adjust the abundances slightly 
in order to coiTect for the highest changes in the depth of the 
edges. Therefore, despite the significance of the different dust 
cross sections, the total opacity is similar to that of pure atomic 
gas and this is the main reason of the only modest improvement 
of the fit with model B compared to model A. It also explains 
why it is hard to use our data to constrain the properties of the 
dust - density, grain size, etc. 




15 20 
wovelength (A) 

Fig. 8. Difference in optical depth for an equivalent hydrogen 
column density of 3.16 x 10^^ m~^, for dust with depletion fac 



tors as given by lWilms et al. I (l200flh . compared to the same ma- 
terial but dissolved into free atoms. The most important spectral 
features are indicated. 



Table 4. Abundances, all relative to the proto-solar values of 
Loddersl(l2003h . 



Element 


Model A 


Model B 


Model C 




pure gas 


gas & dust 


pure gas 




curv. corr. 


curv. corr 


no curv. corr 


H 


= 1 


= 1 


= 1 


N 


1.00 + 0.09 


1.01+0.09 


0.95 + 0.10 





1.012 + 0.014 


1.030 + 0.016 


1.022 + 0.014 


Ne 


1.72 + 0.11 


1.72 + 0.11 


1.83 + 0.11 


Mg 


0.82 + 0.21 


0.85 + 0.21 


0.82 + 0.21 


Fe 


0.76 + 0.07 


0.78 + 0.05 


0.79 + 0.07 



However, the most robust statement we can make about the 
dust is the associated column density. From Table |2] we see 
that the column density of the gas associated with the dust is 
about 70 % of the total column density, and its value of (2.23 + 
0.18)x 10^^ m~^ is in good agreement with the value derived by 



Seward et al.l (l2006h based on modelling of the dust scattering 
halo. Thus, we confirm their finding of a reduced dust column 
towards Crab as compared to a standard gas/dust mixture. 

4.4. ISM abundances 

We have determined the abundances of all elements for which 
we have reliable individual column density measurements. 
These abundances are shown in Table |4] They include all 
phases of the ISM except for the hot ionised gas. We see that 
the abundances do not depend significantly on the inclusion of 
dust in our models. 

Interestingly, the abundances that we derive for most ele- 
ments except for neon are close to the proto-solar abundances 
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of lLoddersI ( l2003h that we used as a reference. Oxygen and ni- 
trogen fully agree with the proto-solar abundances, while the 
magnesium and iron abundances are only slightly lower 

The only exception to this general abundance pa t tern is 
neon. It has an overabundance of ~ 1.7. Antia & Basd (2005) 



and iDrake & Testal (l2005h recently argued that the solar neon 
abundance should be higher. Based on measurements of neon 
to oxygen in stars, Drake & Testa derive a Ne/O number ratio 
of 0.41. Relative to the proto-solar values (Ne/O number ratio 
0.15), this would imply an overabundance of neon compared 
to oxygen of 2.6. Clearly, we find an enhanced neon abun- 
dance, but only 2/3 of the value of Drake & Testa. In fact, our 
Ne/O number ratio of 0.26 ± 0.02 is in excellent agreement 
with recent measuremen ts of the abundances in the Orion neb- 
ula (Ne/0= 0.24 ± 0.05. lEsteban et al.ll200 4') and of B-stars in 
that nebula (Ne/0= 0.26 + 0.07. ICunha et aL.2006.) . 



4.4.1 . The role of carbon 

As the carbon K-edge is longwards of the RGS wavelength 
band, it is difficult to estimate the carbon abundance from our 
measurements. However, as carbon contributes approximately 
10 % to the total opacity below the oxygen K-edge, it is poten- 
tially an important factor to take into account. 

Our best-fit model predicts a total carbon column density 
of 9.2 X 10^' m-2 for our adopted proto-solar abundances. SOO 
found from their optical and UV spectroscopy a column density 
of atomic Ci of 8.5 x 10^' m"^; however the error bars of a 
factor of 2 allow a range between 4- 16x10^' m"^. In addition, 
they found a C ii column of 14x 10^' nT^, but here the error bars 
are a factor of 3, allowing a range between 5 - 50 x 10^' nT^. 
Given these numbers and the fact that a part of the carbon may 
be contained in dust, it is not very likely that the total column 
density of carbon is much below the proto-solar value. 

To estimate the possible effects of a different carbon abun- 
dance, we have determined the best fit parameters in case the 
carbon abundance decreases to 0.7 times proto-solar, a similar 
value as the iron and magnesium abundance towards Crab, 
of our fit then decreases by only 1 .4, the hydrogen column den- 
sity increases by 10^"^ m"^ or 3.2 %, /(5) increases by 0.3 % 
and /(40) decreases by 1.7 %. We will use these numbers in 
our estimate of the systematic uncertainty in the RGS effective 
area. 



5. Discussion of tlie effective area calibration of 
the RGS 

We use our spectral fit to determine the absolute effective area 
of the RGS. This is achieved by determining the power-law cor- 
rection factor f{A) from Eqn. (|5]). As we described in Sect. l3.2l 
we determine this from the values of f{A) at 5 and 40 A. From 
Table |2] we see that we can determine these numbers to an ac- 
curacy of 0.9 and 1.6 %, respectively. The differences of these 
values for model A and B are not very high, at most 2 % at the 
longest wavelengths. 

The differences between model A and B with model C are 
much higher, as model C ignores the intrinsic curvature of the 



underlying continuum. But despite that the curvature correc- 
tion at 40 A (0.3 keV) is 22 % (see Fig.|5]), the enhancement of 
f(40) is only 12 %. This is because also other parameters in the 
fit, in particular the neutral hydrogen column, have changed. 
The lower hydrogen column for model C partially compensates 
for the neglect of the continuum enhancement at low energies. 
Comparing model C with model A, none of the metallic col- 
umn densities neither the abundances differ by more than Icr. 
However, as we have shown in Sect. 12.51 there are strong ar- 
guments for applying the curvature correction. The systematic 
uncertainties on the correction as derived in Appendix B are 
significantly lower than the correction factor itself. We there- 
fore propose to adopt our model B as the model that agrees 
best with all constraints. 

The relatively low statistical uncertainties on f(A) are not 
the most important errors, however There is a range of sys- 
tematic factors that increase the uncertainty on f(A). The most 
obvious one is of course the 10 % uncertainty on the absolute 
flux of the KO I intrinsic continuum that we used (Sect. |2|i. In 
another paper ( Kaastra et al. 2008ah we show how we can re- 
duce that uncertainty by combining the present results from our 
Crab analysis with observations of white dwarfs, for which we 
can determine the absolute X-ray flux with an accuracy of a 
few percent. Other systematic uncertainties cannot be reduced, 
however In order to be able to assess the accuracy of our cali- 
bration, we need to make a quantitative assessment of all these 
systematic effects. 

We estimate the known systematic uncertainties and 
present them in Table |5] We proceed as follows. For each rel- 
evant factor, we estimate its magnitude for a range of charac- 
teristic wavelengths, covering the RGS range between 5-30 A. 
The systematic uncertainties are calculated by comparing the 
nominal fluxes with the derived fluxes when one of the rele- 
vant parameters is changed as described in Appendix B. These 
numbers are given in the table. Then we put each error in one 
of two specific categories. 

For category "c", the systematic uncertainties at different 
wavelengths are more or less correlated. Example: the abso- 
lute flux of Crab is known to only 10 %. Therefore, if the 
true flux of Crab at 5 A differs by 10 % from the value we 
adopted, it will differ by the same percentage for all other wave- 
lengths. Taking this into account, we can assess for these "cor- 
related" errors how they will affect the normalisation and slope 
of f{A) = A{A/ 10)" by "fitting" the estimated systematic un- 
certainties (all having the same sign) directly to a power-law 
shape. 

For category "u", the systematic uncertainties at different 
wavelengths are uncorrected over the full RGS wavelength 
scale. A good example are the cross-dispersion corrections: as 
we use a combination of splines and power-law approxima- 
tions, systematic errors are correlated over ranges between 0.5 
and several A, but there is no relation between the uncertainty 
at 7 A to that at 30 A. In those cases, we estimated AA from 
the typical slope and normalisation that we would obtain if the 
systematic uncertainties at different wavelengths had random 
signs. In general, for this type of uncertainty the direct effect 
on A and a is smaller than for correlated errors (category "c"). 
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A 


E 


a 


b 


c 


d 


e 


f 


g 


h 


i 


j 


k 


1 


m 


A 


keV 




























5 


2.5 


0.1 


0.008 


0.003 


0.000 


0.000 


0.001 


0.001 


0.003 


0.050 


0.000 


0.30 


1.0 


0.0031 


7 


1.8 


0.1 


0.011 


0.005 


0.000 


0.003 


0.002 


0.002 


0.000 


0.010 


0.001 


0.04 


0.022 


0.0025 


10 


1.2 


0.1 


0.015 


0.007 


0.001 


0.005 


0.004 


0.002 


0.004 


0.004 


0.002 


0.01 


0.010 


0.0025 


15 


0.8 


0.1 


0.019 


0.010 


0.001 


0.000 


0.006 


0.010 


0.008 


0.005 


0.002 


0.00 


0.008 


0.0026 


20 


0.6 


0.1 


0.022 


0.011 


0.002 


0.005 


0.007 


0.001 


0.010 


0.015 


0.001 


0.00 


0.008 


0.0045 


25 


0.5 


0.1 


0.024 


0.013 


0.002 


0.005 


0.009 


0.003 


0.012 


0.020 


0.001 


0.00 


0.008 


0.0048 


30 


0.4 


0.1 


0.026 


0.014 


0.003 


0.002 


0.010 


0.001 


0.014 


0.050 


0.000 


0.00 


0.010 


0.0090 




type 


c 


c 


c 


c 


u 


c 


c 


c 


u 


c 


u 


u 


u 




Aa 


0.016 


0.010 


0.006 


0.002 


0.001 


0.005 


0.001 


0.010 


0.010 


0.000 


0.027 


0.006 






AA 


0.100 


0.019 


0.010 


0.001 


0.001 


0.004 


0.005 


0.004 


0.003 


0.002 


0.005 


0.002 





Main power-law: 

" Uncertainty in the KOI continuum model 
Curvature corrections: 

' A 50 % change in the lower energy of the 2-8 keV range used for the normalisation of the Chandra data 
Break size AF = 0.006 at low £ if a break (at 4 keV) is present in the Kuiper et al. data 
0.5 % of the absorption with A'h higher by 7.5 x lO^'* m"^ 
Numerical approximation of fc{E) by ([3} 

Interstellar absorption: 

f Systematic uncertainty in the dust model 

* Neglect of hot gas 

* Interstellar carbon abundance 
Instrumental corrections: 

' Cross-dispersion correction 

^ Pile-up correction for 2 CCD mode; adopts 5 % uncertainty in a (la) 

* Redistribution errors, estimated from difference fit residuals RGSl and RGS2 
' Errors in the Chebychev polynomials used for the effective area correction 

Typical combined RGS1/RGS2 statistical uncertainty for a full 1 A wide bin 



More details on the individual systematic uncertainties are 
given in Appendix B. 

For completeness, we also list in Table|5]the statistical un- 
certainties on the Crab spectrum. As can be seen, these are low 
compared to the systematic uncertainties. 

The main systematic uncertainty in the slope of our power- 
law coiTection (+0.027) originates from the redistribution er- 
rors (column (k) of Table |5]l. These uncertainties are mostly 
important below 10 A. In principle, it should be possible to re- 
duce these errors further, by a carefull analysis of a large sam- 
ple of sources with different intrinsic spectral shape in the hard 
band. However, such an effort is beyond the scope of this paper. 
If we add up all other systematic uncertainties Aa on the slope 
in quadrature, except for the redistribution errors, the combined 
uncertainty is 0.020; the total rm.s. systematic eiTor is 0.034. 

The main uncertainty on the absolute flux is - apart from 
the uncertainty in the KOI model listed in column (a) - the 
uncertainty in the lower energy used for the scaling of the 
Chandra spectra on the KOI model (column (b). 



In another paper dKaastra et alj|2008a[) we show how this 
uncertainty in flux can be reduced to about 3 % by combining 
our results for Crab with white dwarf spectra. However, none 
of that analysis affects our present results on the Crab nebula. 



6. Conclusions 

We have studied the interstellar absorption towards Crab and 
determined the effective area of the RGS. 

The interstellar line of sight towards Crab appears to be free 
of molecules and has a relatively low contribution of hot gas. 
It contains only about 70 % of the amount of dust expected 
from its absorption column, a conclusio n that is in agreemen t 
with studies of the dust scattering halo bv lSeward et al. I (l2006h . 
From the centroids of the dominant O-K and Fe-L absorp- 
tion features, we conclude that gas and dust contribute roughly 
equally to these features, and that ferric iron is more abundant 
in the dust than ferrous iron. 

We have determined the interstellar abundances with very 
high precision, the oxygen abundance with a relative accuracy 
even better than 2 %. While oxygen and nitrogen are fully 
compatible with proto-solar abundances, we find slightly lower 
(~ 80 % proto-solar) abundances for iron and magnesium, al- 
though for the last element we cannot exclude proto-solar abun- 
dances. On the other hand, neon is over-abundant by a factor of 
1.7. 

We can also put tight limits on the amount of singly ionised 
ions, thanks to our high spectral resolution. For instance, less 
than 2 % of all oxygen is in the form of O ii. 

Concerning the calibration of RGS, we have accurately de- 
termined the effective area of this instrument. The relative accu- 
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racy of the effective area over the 7-30 A band is equivalent to 
a power-law index uncertainty of about 0.03. The absolute ef- 
fective area from our analysis is known to the same accuracy as 
the adopted continuum model s pectrum of Crab, namely 10 %. 
In a forthcoming paper (iKaastra et al.l2008ai) we will show how 
that last uncertainty can be reduced to a few percent. 
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Appendix A: Spectral extraction and response 
generation 

A.1. Spectral extraction 

Spectra and response matrices were created using the stan- 
dard XMM-Newton data analysis system (SAS), version 7. The 
standard RGS pipeline 'rgsproc' was used, but with a few adap- 
tations explained below. 

For point sources, the preferred method is to extract spectra 
that contain 95 % of all photons in the cross-dispersion direc- 
tion, and 95 % of all photons in the pulseheight-wavelength 
(PI, "banana") regions. The response matrix then corrects for 
the missing flux fraction in both subspaces. 
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In our case, the situation is different. Crab is an extended 
source, with substantial emission wings containing ~10 % of 
all flux in the ROSAT band. These wings a re visi ble up to 
half a degree from Crab (Predehl & Schmitt I l995h . Crab is 
therefore larger than the width of the RGS detector in the 
cross-dispersion direction, and also photons are spread over a 
larger region in the dispersion direction than for a point source. 
Moreover, Crab is bright and for some parts of the source pile- 
up of events cannot be neglected. 

For these reasons, we extracted source spectra over the full 
width of the RGS detector (5.0 arcmin), and with a PI selection 
that contains 98 % of the flux for each wavelength. Selecting 
100 % P/ is not possible because then spectral order separation 
becomes impossible and at low PI values the selected range 
would fall below the lower threshold. 

To prevent the broader PI selection to select background 
events from the CCD system peak, the PI selection was limited 
to PI values higher than 200 eV, well above the average system 
peak of around 100 eV. This made the PI selection narrower 
than 98 % for wavelengths longer than 24.8 A. This limitation 
of PI selection for a particular wavelength range is properly 
taken into account by the SAS processing system. 

The response matrices that we used, however, were created 
for a selection of 95 % in PI space and 95 % in cross-dispersion 
space. These are the recommended values for point sources. 
We use those matrices because they are the best calibrated (see 
Sect. IA.2I ). but we apply correction factors in the form of an 
arf file to this standard response matrix. Note that in general 
the standard SAS software produces only rmf matrices where 
the effective area has already been included. As our corrections 
go beyond the standard SAS, we have chosen here for the arf 
option. 

The correction factors entering the arf are discussed in the 
next sections. All factors are multiplied to get the total area 
correction. They are: 

1. 0.98/0.95, to correct for the 95 % PI selection in the stan- 
dard response matrix as compared to our actual 98 % ex- 
traction. 

2. 1.00/0.95, to correct for the 95 % cross-dispersion selec- 
tion in the standard response matrix as compared to our full 
detector (100 %) extraction. 

3. Dead time can be ignored. See Sect. |A.3l 

4. Pile-up correction /p(/t). See Sect. |A.4| 

5. Vignetting of telescope and gratings. See Sect. |A.5l 

6. Loss of photons in the cross-dispersion direction due to the 
extended size of Crab /d(/l). See Sect. lA.6l 

Finally, we need to correct for the extent of the source in 
the dispersion direction. This causes some spectral resolution 
degradation, but does not affect the effective area. This is dis- 
cussed in Sect. IA.7I 



A.2. Standard response matrix 

The standard response matrix contains an empirical correction 
factor for the nominal effective area for point sources (95 % 
PI selection, 95 % cross-dispersion selection). This empirical 



correction is based on all RGS spectra from Mrk 421, used as 
calibration source. We assume Mrk 421 to have a pure power- 
law spectrum of arbitrary power-law index and normalisation, 
which may change over time. We fit a power-law in the 10- 
25 A wavelength range. The remaining part of the RGS spec- 
trum is made to fit this spectrum using a set of low frequency 
Chebychev correction polynomials (up to order 12). A time- 
dependant correction is added assuming an increasing Carbon 
deposit on the RGS CCD's, with parameters based on the con- 
stant sources RX J1856. 5-3754 and the Vela pulsar. 

The standard SAS calibration now contains the Chebychev 
correction polynomials and carbon absorption model including 
the absolute effective area calibration based on the work de- 
scribed in this article. This latter correction was of course not 
used for this analysis, since this is a factor to be determined. 

A.3. Dead time effects 

For sources producing high count rates in the RGS such as the 
Crab nebula, many events are lost due to onboard processing 
constraints. When too many events are queued for on board 
processing or telemetry download, the processor will discard a 
number of CCD readouts (frames) to free its buffer. This will 
influence dead time, but will have no other consequences for 
the instrument response. The data reduction SAS software will 
take the exact number of lost frames into account to compute 
the proper exposure time. Net result will only be a higher sta- 
tistical noise than is expected on the basis of observing time 
only. 

Since we selected 100 % of the CCD width, events recorded 
during CCD frame transfer (out-of-time events) would still fall 
within the spatial selection window and hence do not contribute 
to any dead time. 

A.4. Pile-up corrections 

The brightest parts of the Crab nebula suffer from pile-up in the 
spectrum. This is evident from a reduced count rate in observa- 
tions with 8 CCDs readout as compared to single CCD readout 
mode. For a given CCD, the net exposure time per frame in k 
CCD readout mode is simply k times higher than the default 
integration time of 0.57 s. 

It is difficult to make an exact model for the pile-up. In 
general, any photon hitting the detector may give a signal ei- 
ther in a single pixel or in more than one pixel, with a range of 
many possible patterns. Within the same (time) frame, a second 
photon may hit the same or a neighbouring pixel, with its own 
pattern. The combined number of possible patterns is large, and 
the electronics, on-board and ground software treat each com- 
bination differently, also depending on the energy of both pho- 
tons. 

Therefore we create an empirical pile-up model. We de- 
fine fi to be the expected number of counts per frame per pixel 
in single CCD readout mode (in the absence of pile-up). This 
number is obviously a function of the position on the detector. 
Even for sources as bright as Crab, the count rates are mod- 
erate. The maximum value of yu for the Crab nebula is less 
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than 0.004. Therefore, already double events are rare and we 
can fully ignore pixels with three or more events. We consider 
only the effect of pile-up on the first-order spectrum. In general, 
when a genuine 1st order photon hits a pixel, and within the 
same frame a second event (first, second order or even back- 
ground) hits the same or a neighbouring pixel, several things 
may happen but all lead to a rejection of the event as a first or- 
der photon: it is either rejected completely, or seen as a second 
order event with twice the energy or alternatively as a back- 
ground event. Accounting for pile-up, yUi p, the expected num- 
ber of first-order events in a given pixel can be written as: 

fJ^\,p = - ak/j.) (A.l) 

where and /u are the expected first order and total counts/ 
per pixel per frame without pile-up, k is the number of CCD's 
read-out and a is a constant that needs to be determined. 
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Fig. A.l. Ratio of first order count rates for the Crab nebula for 
8 CCD mode as compared to single CCD mode, as a function 
of fi, the total number of counts per frame per pixel. We only 
show the results for CCD 5-8 of RGS2. The results for the 
other CCDs are similar The curve indicates our best-fit model 
to Eqn. (lA.ll i. 



Fitting the pile-up correction in Eqn. iA.li not to (1 - akfi) but 
to {b - akfj), we find that b is typically a few percent lower 
than 1 . However, as the fit is statistically driven by the brighter 
parts with high count rates (see the much lower error bars for 
high values of //), the effect on the derived pile-up correction at 
high count rates is small. These brighter parts also dominate the 
total count rate of Crab. We estimate that the net effect on the 
integrated spectrum is equivalent to an enhanced uncertainty of 
~ 0.45 in a, i.e. about 3 times the nominal statistical uncertainty 
in that quantity. 

The value for a can also be used for observations of other 
targets, but it should be remembered that in our case it is de- 
fined for the 98 % selection in energy space that we needed to 
use because of the extended nature of Crab. 

The value of a = 6.45 is used to derive a wavelength- 
dependent correction factor for pileup fp{A), that we for- 
mally include in the effective area correction ("arf ). It is con- 
structed by averaging 1 - akfj. over the cross-dispersion direc- 
tion, weighted by the observed first order count rate image cor- 
rected for pile-up. 
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Fig. A.2. Pile-up correction for the Crab nebula observed with 
RGS2 in 8 CCD mode. Most of the spikes of lower pile-up are 
caused by CCD gaps and the central columns of CCDs. 



For each pixel and all CCDs we have estimated /i by aver- 
aging all counts over time and 7x7 pixels around the relevant 
pixel (in order to reduce the noise). Then we grouped the pixels 
together in groups with the same range of fi values, and deter- 
mined the total number of first order counts for the group in 
each readout mode. The number of first order counts for differ- 
ent k, as a function of fi, were fitted to Eqn. (lA.ll ). An example 
is shown in Fig. lA.ll We determined the weighted average of a 
from all CCDs. These values are 8.17 + 1.17, 6.46 + 0.30 and 
6.41 + 0.17, from a comparison of A: = 1 with k - 2, k = 4 
and ^ = 8, respectively. Combining these numbers, our best 
estimate of a is 6.45 + 0.15. 

There is no significant difference between the value of a for 
RGSl and RGS2. As can be seen from Fig. I A. II at low count 
rates there is a small, unexplained, deficit of count rate ratio. 



In Fig. lA.2l we show the arf correction for the Crab nebula. 
In 8 CCD mode, the maximum correction occurs at 1 keV and 
amounts to 14 %. 

A.5. Vignetting 

The effective area of the RGS depends on the position of the 
source in the field of view. The standard effective area is cal- 
culated for a point source. For off-axis photons (either in the 
dispersion or cross-dispersion direction), the effective area dif- 
fers due to vignetting of the telescope and obscuration and vi- 
gnetting by the gratings. Crab is an extended source. For each 
energy band we have calculated the effective vignetting by mul- 
tiplying a MOS image (see Sect. IA.7b with the vignetting fac- 
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tor, including both telescope vignetting and grating obscura- 
tion. The average vignetting factor is 0.983, and to high accu- 
racy it does not depend on energy within the RGS wavelength 
band. There is a weak increase of 0.2 % in telescope vignetting 
from 7-25 A caused by the enhanced size of Crab at longer 
wavelengths, and therefore more photons are off-axis. But this 
increase is compensated for by slightly lower grating obscura- 
tion at longer wavelengths. Therefore we multiply our effective 
area by this constant factor of 0.983. 

A.6. Cross-dispersion correction 

As Crab is an extended source and the RGS detectors have a 
limited width of about 5' in the cross-dispersion direction, a 
fraction of all photons at each wavelength will not hit the detec- 
tor. There are three causes for this loss: the extended emission 
of the nebula, dust scattering and additional scattering by the 
gratings of the RGS. In order to be able to use spectral mod- 
els for the full Crab nebula, we need to account for this loss of 
events. We have implemented this as follows. 

RGS2 CCDS 
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Fig. A.3. Cross dispersion profile for CCD6 of RGS2 for our 5 
observations (central and two off-axis pointings at each side). 
The data have been integrated over the full CCD in the dis- 
persion direction, and over 8 pixels in the cross-dispersion di- 
rection. The best-fit power-law models, fitted individually to 
chipy< and chipy> 100, are shown by the thick solid lines. 
Note that we aligned all off-axis pointings to the on-axis (chipy) 
coordinate in the cross-dispersion direction. The different sym- 
bols correspond to the different observations. 



For our spectral fitting, we only use the on-source pointing, 
but we use our off-axis pointings in the cross-dispersion direc- 
tion to estimate for each wavelength interval the fraction of the 
total flux that is lost in the central pointing. Unfortunately, the 
extent of the Crab nebula is even larger than our maximum off- 
axis angle (9.5')- Therefore we have to extrapolate the flux out- 
side this maximum off-axis angle. In general, power-laws are 
good approximations for the dust scattering halo's at a given 



energy jPredehl & Schmittlll995l:IPredehl & Kloselll996h . Due 

to the asymmetry of the nebula and the grating scattering, we 
determine separate correction factors for both sides of the neb- 
ula. For the longest wavelengths, background becomes impor- 
tant and we have modelled this by subtracting an observation 
of the Lockman hole. 

A typical example of our fits is shown in Fig. lA.3l It is seen 
that the individual pointings do not fully coincide at their low- 
est and highest chipy positions. This is due to a combination 
of the outwardly decreasing flux of the nebula with the slightly 
asymmetric scattering of the gratings. However, the total cor- 
rection factor is hardly affected by this slight lack of overlap, 
as we use two full off-axis observations on each side and our fit 
averages out these boundary effects. 




15 20 
Wavelength (A) 

Fig. A.4. Enclosed flux fraction for the on-axis pointing on the 
Crab as a function of wavelength. Circles are RGSl data, and 
triangles RGS2 data, shifted by -1-0.02 in enclosed fraction. The 
solid curve shows our finally adopted smoothed approxima- 
tion to this correction factor; the dotted line shows the original 
spline fit to the data as described in the text. 



We determined the cross-dispersion loss correction factor 
/d (defined as the fraction of the Crab flux falling on the RGS 
detector for the central pointings) with the highest possible res- 
olution. In practice this is achieved by dividing each CCD into 
8 equal wavelength bins. However, in particular for the CCDs 
with lower count rates, we had to go to 4 or even 2 bins per 
CCD in order to obtain accurate values. We show the mea- 
sured correction factors in Fig. IA.4I The following should 
be noted. 

First, for RGS2 /d is on average 0.02 lower than for RGSl, 
caused by stronger scattering in RGS2. For this reason we have 
shifted in Fig. IA.4I the RGS2 data points by H-0.02, in order 
to eliminate in first order the differences between both RGS 
detectors. 

Secondly, the curve shows fine structure near some impor- 
tant edges, such as the O-K edge (23 A), Fe-L (17 A) and Mg K 
(9 A) edges. This structure is real and is caused by the dust scat- 
tering, which is wavelength dependent and shows structure at 
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these edges (see for example lDrainell2003l and lCostantini et al 
2005b . Unfortunately we have insufficient statistics to increase 
the wavelength resolution of /j near these edges. It is also 
hard to model the fine structure theoretically, as it depends 
on the chemical and size distribution of the dust grains, the 
spatial distribution of the underlying synchrotron emission, 
and the complex redistribution of photons in the combined 
dispersion/cross-dispersion plane of RGS. 

Finally, the accuracy of /d decreases rapidly towards longer 
wavelengths, and beyond 25 A it is impossible to determine 
it at all. This is caused by an accumulation of several ef- 
fects, namely (1) the decreasing flux of the Crab nebula to- 
wards longer wavelengths due to interstellar absorption, (2) 
the stronger dust scattering at long wavelengths, making our 
power-law fits to the spatial profile flatter and hence the correc- 
tion factors higher and more uncertain, and (3) at the low flux 
levels of the off-axis pointings at long wavelengths, systematic 
background uncertainties emerge. 



We approximate the data shown in Fig. IA.4l as follows. As 
a first approximation, we fit a power-law to the scattered flux 
fraction, an assumption that is often made and is based on the 
Rayleigh-Gans approximation but which however may not be 
fully appropriate at the longest wavelengths (see for example 
Smith & Dweklll998l) . The power-law is given by /d - 1 - 
0.00338/i'^^^^ with A the wavelength in A. A part of this tail is 
still visible in our model beyond 27 A, and lacking other data, 
we use this approximation for all A > 21 A. For the remaining 
part, we fitted a spline to the residuals of the power-law fit, and 
added that spline to our correction factor Around 23 and 24 A, 
this spline shows a strong wiggle caused by the two RGS 1 data 
points at those wavelengths (see the dotted curve in Fig. lA.4b . 
Such a wiggle is physically not expected and it is perhaps due 
to a little overshooting of our fit near the end points. It also 
appears that after spectral fitting the predicted flux in the 1- 
2 A wide band just longwards of the O-K edge was too high 
by up to 4 %. Therefore we have fitted a spline to the spectral 
fit residuals of our preliminary fit in the 23-26 A range, and 
included that spline in our final approximation for /d, the solid 
curve of Fig. lA.4l As our spectral fits show (see Sect. 13.31 ) this 
gives reasonable results. 



A.7. Broadening in the dispersion direction 

The Crab nebula is an extended source. The first order spec- 
trum of any part of the remnant, at projected distance AO (in 
arcmin) along the RGS dispersion axis, will be shifted in wave- 
length by AA (A) = 0.138Aft This holds for all wavelengths. 
Therefore, the intensity profile I{6) of Crab, integrated over 
the cross-dispersion direction for the region covered by the 
RGS detector, translates directly into a wavelength broaden- 
ing profile with the same shape. In our spectral fitting, after 
deriving the source spectrum S (A) we convolve this with the 
broadening kernel I{AA) before folding with the response ma- 
trix. Technically, this is achieved by using the Ipro convolution 
model of the SPEX package, which has been used successfully 
for observations with RGS of sever al extended X-ray so urces 
like clusters of galaxies, for example iKaastra et alj (1200 lb . 



An underlying assumption of this broadening model is that 
the spatial profile is independent of wavelength. Our discus- 
sion of the Chandra imaging data of M04 (Sect. 12.5.21 ) shows 
that this is strictly speaking not true for the power-law contin- 
uum. However, for the continuum the effect of the broadening 
is not important, as it varies on much larger wavelength scales 
than the broadening kernel, and M04 also showed that there 
is little spatial variation in the interstellar absorption column. 
Therefore our approximation can be used safely. 




Fig. A.5. Broadening kernel I(AA) to be used to model the spa- 
tial extent of Crab in the dispersion direction. The kernel is 
centred around the nominal position of the pulsar The average 
wavelength shift of the kernel is located at AA = -1-0.014 A and 
is indicated with the dashed line. 



We have constructed the broadening kernel I(AA) using an 
XMM-Newton observation taken in March 2002 (orbit 4 11 , ob- 
sid 0153750501). That observation had the same pointing cen- 
tre and roll angle as our present observation. MOS 1 observed 
the source in full window mode. The central CCD of MOSl 
suffered from severe loss of events due to telemetry limitations, 
which are more important than pile-up for this source. These 
telemetry losses do not affect the spectral shape. The remaining 
fraction of events for the central MOS CCD (3.64 %) could be 
estimated by constructing the projected intensity profiles along 
each axis and matching that to the projected profiles of the outer 
CCDs, that do not suffer from telemetry saturation or pile-up. 
Only the pulsar region suffers from severe pile-up (Fig. IA.6I 1. 
but in the RGS band the pulsar contribution is modest (typically 
3-5 %). This MOSl image in the RGS band (0.35-2 keV) was 
then integrated over the cross-dispersion of RGS, yielding our 
kernel I(AA). The kernel is shown in Fig. lA.5l 

The full-width at half maximum of the kernel is 0.22 A, 
which therefore is our effective spectral resolution. Our spectra 
were extracted with respect to the position of the pulsar. The 
nebula shows a small asymmetry with respect to that position: 
the average centre of Crab is shifted by 0.014 A. 
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Fig. A.6. MOS 1 image of Crab. Shown is a part of the central 
CCD image of observation 0153750501, on a linear intensity 
scale. North is up. RGS spectra are extracted between the two 
parallel lines (separated by 5 arcmin). The pointing is centred 
on the pulsar, which shows pile-up in this image. 
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Fig. B.l. Relative residu als of the BeppoSA X measurements of 
the Crab nebula used by lKuiper et al.l (1200 ih (Fig.|2]i compared 



to the best fit model (F = 2.147 ± 0.002 and = 9.31 ± 0.05). 
For clarity, error bars have been omitted, but the solid lines in- 
dicate the + 1 o" statjsticaj error range. Circles: MFCS; triangles: 
HPGSPC; diamonds: PDS. 



A.8. Spectra used for the spectral fitting 

Finally we need to decide which data sets to use. We have 
spectra in 1, 2, 4 and 8 CCD readout mode. The pile-up cor- 
rections increase with the number of CCDs that are readout 
simultaneously. As a result, the fits for 8 CCD mode may have 
slightly higher systematic uncertainties than the other spectra. 
Moreover, the exposure times differ for the readout modes. We 
have made fits for all combinations 1, 1h-2, 1-h2-h4 and 
1 H- 2 H- 4 -H 8 CCD readout mode. As a typical example, the 
best values for the hydrogen column in units of lO^'' m"^ are 
3 . 1 74 ± 0.036, 3 . 1 75 + 0.022, 3 . 1 84 ± 0.020, and 3 . 1 60 + 0.020 
respectively (for model A, without dust). Including the 4 CCD 
mode data does not improve the statistics significantly, while 
including the 8 CCD mode data puts the best fit parameter near 
its Icr limit, indicating that the value for this mode alone devi- 
ates more. For these reasons, we here consider only the com- 
bined dataset of single and double CCD readout. 

Appendix B: Breakdown of systematic 
uncertainties 

In this appendix we discuss in more detail the various system- 
atic uncertainties on our derived effective area solution. 



B.1. Systematic uncertainties on the adopted 
spectrum of Crab 

B.l .1 . Systematic uncertainties on tine absolute flux 
level of Crab 

We adopt here the typical uncertainty on the KOI continuum 
level of 10 %. In addition, there is a small uncertainty related 
to the statistical error of the slope in the KOI data (Sect. |2|i. 



However, there are also systematic uncertainties in the shape 
of the KOI continuum. These are discussed in more detail in 
the next subsection. We list these uncertainties as column (a) in 
Tabled 

B.1 .2. Systematic uncertainties on the shape of the 
Crab spectrum 

The KOI model that we use for the high energy spectrum of 
Crab is based on BeppoSAX data. The statistical uncertainties 
can be determined well (Sect. |2]i, but it is harder to determine 
the systematic uncertainties. 

For the RGS band, t he MFCS dat a are m ost important so 
we consider those first. iBoella et"al 



1997 ') describe in de- 



tail the calibration of the MFCS, based on ground measure- 
ments and careful modelling of the components of the in- 
strument. Unfortunately, no estimates of systematic u ncertain- 



ties a re given. However, the BeppoSAX cookbook (Fiore et al 
1999!) gives some clues. 



First, when the total (nebula plus pulsar) Crab data for 
MFCS, HPGSPC and PDS are fitted individually, the photon 
indices derived from the individual instruments agree to within 
±0.016. 

Secondly, the BeppoSAX cookbook indicates that the aver- 
age and dispersion of the PDS to MFCS normalisation factor 
for five different sources excluding Crab, is 0.86 and 0.03, re- 
spectively (compare to the value of 0.87 obtained by KOI for 
Crab). Adopting this scatter of 0.03 as a typical slope uncer- 
tainty over the effective centres of the MFCS and PDS bands 
at 4.5 and 45 keV, respectively, we find a corresponding maxi- 
mum systematic photon index uncertainty of 0.015. 

Third, we show in Fig. lB.ll in more detail the residuals of 
the fit presented in Fig. |2] These residuals, corresponding to 
calibration uncertainties within the instruments, are typically 
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lower than 5 % and therefore contribute to systematic uncer- 
tainties on the photon index over the 2-100 keV band of at 
most +0.012 log 1.05/ log 50). 

Finally, the latest calibration of INTEGRAL, covering also 
a broad energy range just like the BeppoSAX data used by 
KOI, gives similar results for the hard X-ray part of the spec- 
trum (IJou rdain et al. 2008). Down from 3 keV up to the spec- 
tral break at 100 keV, they fit a single power-law with F = 
2.105 + 0.030 to the full Crab (nebula plus pulsar) spectrum, 
using exactly the same Galactic foreground absorption model 
as KOI. Over the 2-8 keV band, the KOI model (curved pulsar 
plus nebula spectrum) has an effective photon index of 2.109, 
only 0.004 higher than the INTEGRAL value. Moreover, the 
INTEGRAL flux in this 2-8 keV band is 6 % higher than the 
KOI flux in the same band, well within the quoted 10 % uncer- 
tainty. Thus, if we would have adopted the INTEGRAL model 
spectrum instead, the main change would be an overall ~ 6 % 
higher flux hence a ~ 6 % lower eff'ective area for RGS, but no 
significant slope diff'erence. 

The above estimates can be summarised in a maximum sys- 
tematic slope uncertainty of about 0.016. 

B.2. Systematic uncertainties on tine curvature 
correction 

The curvature coiTection based on Chandra imaging of the neb- 
ula as derived in Sect. l2.5l has its limitations, leading to system- 
atic uncertainties in its size. We assess these uncertainties here. 

B.2.1 . Energy range normalisation Cliandra continuum 

The choice of the 2-8 keV band to normalise the continuum af- 
fects the correction slightly. Changing the lower limit by 50 % 
to 1.5 or 3.0 keV modifies fdE) typically by less than 2.6 %. 
We list this as column (b) in Table|5] 

B.2.2. Local high-energy cut-offs 

Some of the local power-law fits of M04 may be affected by 
unnoted high-energy cut-offs. If this happens, a power-law fit 
gives a too soft spectrum, resulting in an over-estimate of fdE) 
at low energies. This effect is hard to quantify, but the fact that 
the photon index of F = 2.10 derived from a direct fit to the 
total ACIS-S spectrum (M04) agrees well with the spatially in- 
tegrated spectrum indicates that it is not important. Also, if a 
significant fraction of the pixels would have a high-energy cut- 
off, a clear break in the BeppoSAX data used by KOI should be 
present, which is not the case. The maximum change in pho- 
ton index for an assumed spectral break near 4 keV in these 
BeppoSAX data is less than 0.006. This corresponds to a max- 
imum error in fc of at most 1.4 % (column (c)). 

B.2.3. Pile-up corrections Chandra 

M04 needed to make pile-up corrections; they estimate a sys- 
tematic uncertainty of 0.05 in their photon indices. As spectral 
curvature corresponds to differential photon index, the effect is 



small (less than 0.03 %) and we shall ignore it here. Another 
potential effect might be a systematic error in the pile-up cor- 
rection of M04 that would correlate with intensity and photon 
index. The pile-up correction for the photon index F of M04 is 
always less than 0.17 forF =1.5 and less than 0.21 for F = 2.5. 
We may assume very conservatively that the maximum system- 
atic error is half of this difference, i.e. 0.02. If we increase and 
stretch all local photon indices F with AF(F) - 0.02(F-2. 1) and 
re-evaluate the curvature coiTection, we find differences lower 
than 0.6 % for A < 30 A. As for most wavelength the cor- 
rection is much lower and the maximum of 0.02 is likely an 
over-estimate, we can safely ignore this effect here. 

B.2.4. Different continuum model for different 
instruments 

In the 1-8 keV band, there is a difference of up to 10 % between 
the model flux of KOI and M04, and a corresponding difference 
in photon index of 2.147 - 2.10 = 0.047. Such differences are 
typical for spectra taken by different instruments, and generally 
they coiTespond to calibration uncertainties. As above, we do 
not expect this difference to produce artificial spectral curva- 
ture. 

B.2.5. Interstellar absorption column density 

There is a big difference between the interstellar absorption A'h 
used in both papers: KOI have A^h = 3.61 x 10^^ m"^, while 
M04 find 3.20 x 10^^ m~^. In the absorbed model continuum 
near the interstellar oxygen edge at 0.53 keV, this causes a 
30/50 % higher flux below/above the edge in the M04 model 
as compared to the KOI model. Apparently, the difference is 
caused by a mismatch of the same order of magnitude between 
the effective areas of BeppoSAX/LECS and Chandra/ACIS-S 
near the oxygen edge. Nevertheless, this does not affect the 
high-energy continuum, and it is also not too important for the 
determination of the spectral curvature, as that is caused by true 
spatial variations of the photon index. 

B.2.6. Biased absorption due to region selection 

The value of A^h derived by M04 may be slightly biased, be- 
cause in order to avoid pile-up, it was determined in the low 
surface brightness regions. These regions are mostly found in 
the outer parts and have softer intrinsic spectra. More impor- 
tantly, the softest X-ray flux in these regions may be enhanced 
by up to ~ 10 % due to dust scattering from the brighter cen- 
tral parts. This may cause a slight under-estimate of the derived 
interstellar column, as dust scattering is ignored. On the other 
hand, the same regions may loose soft X-ray flux to regions 
even further out, thereby reducing the effect. It is therefore hard 
to quantify this effect, and we will ignore it. 

B.2.7. Spatial variations of Nu 

M04 assumed that all parts of the nebula have the same inter- 
stellar column density. This is a reasonable assumption, as they 



J.S. Kaastra et al.: X-ray spectroscopy of the Crab nebula 



23 



show that less than 0.5 % of all pixels has excess absorption up 
to 7.5 X lO^'* m"^ above the average value. The excess is caused 
by filaments visible in optical line emission. The effect of this 
on the integrated spectrum is less than 0.3 % near the oxygen 
edge. We list this as column (d) in Table|5] 

M04 also estimate that the scatter in A^h due to spatial varia- 
tions of the interstellar column density is less than 2 x 10^"* m"^. 
For the integrated nebular spectrum, most of these variations 
will average out; ignoring systematic effects, the uncertainty 
in the mean absorption for the 2074 pixels considered by 
M04 would be ~ V2074 lower than the scatter, hence about 
4 X 10^^ m"^, which is undetectable. 



B.4.3. Redistribution errors 

We have calibrated the effective area by multiplying the nom- 
inal effective area by an empirical fudge factor. However, at 
short wavelengths important uncertainty is introduced by wide 
angle scattering of photons by the gratings of RGS. This can- 
not be described by a simple arf-correction, as it depends on the 
spectral shape of the source. An estimate of the corresponding 
uncertainty in the effective area is obtained by comparing the 
fit residuals for RGS 1 and RGS2 separately. We take the differ- 
ence in residuals, averaged over 1 A wide bins, as an estimate 
for the systematic uncertainty associated to this effect. The un- 
certainties are shown in Table|5] column (k). 



B.3. Systematic uncertainties in tlie absorption model b.4.4. Errors in tine Cliebycliev polynomial corrections 



B.3.1. The role of dust 

There are small differences between the continuum model de- 
rived from model A (without dust) and model B (with dust). 
In principle, model B is our preferred model but it still may 
have its systematic uncertainties, related to issues as the exact 
density, grain chem ical composition and size distribution. As 
Seward et al. I (l2006h typically found differences up to 30 % in 
the amount of dust towards Crab depending on their adopted 
dust model, we will, quite arbitrarily, take here half of the dif- 
ference between model A and B as the systematic uncertainty 
associated to the dust model. This is listed in column (f) of 
Table m 



B.3.2. Hot gas 

The maximum amount of hot gas in the line of sight towards 
Crab is A^H ^ 5 X 10-^ m"- (Sect. 14.3. li We take the transmis- 
sion of gas with such a column density, proto-solar abundances 
and a temperature of 0.2 keV as our estimate for the associated 
systematic uncertainty introduced by ignoring this component 
(column (g) of Table|5]l. 



Our initial effective area was based on observed spectra of 
Mrk 421; mismatches between the observed spectra and the 
adopted power-law continuum of these sources were corrected 
using a Chebychev polynomial expansion of a correction fac- 
tor for the effective area of RGS (see Sect. IA.2I ) . We estimate 
the uncertainty associated to this by analysing the fit residu- 
als of another source, PKS 2155-304, orbit 362. The resulting 
uncertainties are shown in column (1) of Table |5] 

B.4.5. Errors in adopted spatial profile in the 
dispersion direction 

The Crab spectrum is energy dependent. We have used a MOS 1 
image in the 0.35 - 2.0 keV range to predict the spatial distri- 
bution of the emission in the dispersion direction of RGS. We 
have chosen a number of narrow band energy intervals in the 
MOS image to see how much that projected profile deviates 
from the mean 0.35 - 2.0 keV profile. These intervals range 
from 8 to 22 A. We find shifts lower than 1 mA in the centroid 
(to be compared to the FWHM of the kernel of 220 mA), and 
no deviations higher than 2 % in the FWHM of the profile. This 
justifies our choice of a single broadening kernel. 



B.3.3. The carbon abundance 



We follow here Sect. I4.4.T] and adopt a 30 % systematic uncer- 
tainty in the interstellar carbon abundance. See further column 
(h) of Table |5] 



B.4. systematic uncertainties related to the instrument 

B.4.1 . The cross-dispersion correction 

Our estimate is based on the residuals shown in Fig. |A.4| and is 
shown in column (i) of Table |5] 



List of Objects 

'Crab' on page[l] 
'Cyg X-r on page [13 
'4U 1820-303' on pageE 
'Mrk 421' on page[l7 



'RX J1856.5-3754' on 
'Vela pulsar' on page 
'PKS 2155-304' on page|23| 
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B.4.2. Pile-up correction 

Our estimate is based on Fig. |A.2| adopting a 5 % (2cr) uncer- 
tainty in the pile-up scale parameter a. The uncertainties are 
shown in Table|5] column (j). 



